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Nomenclature 


3- vectors are denoted by bold face quantities such as p. 

4- vectors are denoted by non-bold face quantities such as p. 
cm is the abbreviation for the center of momentum frame. 

The lab frame is always identical to the target frame, which is typically the spacecraft frame.. 

* refers to either the cm or projectile frames. 

Subscript c refers to the cm frame. 

Subscript p refers to the projectile frame. 

Subscript l or lab refers to the lab frame. 

If equations do not have & *, c, p, or l subscript, then it means the equation is true in any 
frame. 

E* or p* or 0* refer to the energy or momentum or angle of an arbitrary particle in the * 
frame. 

E c or p c or 9 C refer to the energy or momentum or angle of an arbitrary particle in the cm frame. 

E p or p p or 9 P refer to the energy or momentum or angle of an arbitrary particle in the projectile 
frame. 

Ei or p^ or 9[ refer to the energy or momentum or angle of an arbitrary particle in the lab 
frame. 

In general, x* refers to the value of the x quantity of an arbitrary particle in the * frame, 
and similarly for the cm and projectile frames. 

Ej* or pj* refers to the energy or momentum of particle j in the * frame, and similarly for 
the cm and projectile frames. 

In general, Xj * refers to the value of the x quantity of particle j in the * frame, and simi- 
larly for the cm and projectile frames. 

9jk* refers to the angle of the momentum of particle j with respect to the momentum of 
particle k. The angle is measured in the * frame, and similarly for the cm and projectile frames. 
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Often, just Oj or 9ji or 6j c will be written. This is used when we are referring to the angle 
of particle j but have not yet needed to decide what the angle of particle j is with respect to. 

Consider 1 + 2 — ► 3 + 4. In the cm frame, Pi c + P2c = 0 = P3 C + P4 C implying |pi c | = |p2c| = |Pic| 
and | P3 C | = |p4c| = |Pf c |- -E'icj E fc or p j c , p f c refer to the energies or momenta of the initial (i.e. 
particles 1 or 2) or final (i.e. particles 3 or 4) particles in the cm frame. 

/3* or 7* refer to the speed (in units of c) or relativistic 7 factor of an arbitrary particle in 
the * frame. They are related via 7 2 = 1 _ 1 o 2 , and similarly for the cm and projectile frames. 

/3j> or 7j* refer to the speed or relativistic 7 factor of particle j in the * frame, and simi- 
larly for the cm and projectile frames. 

/ 3*i or y*z refer to the speed or relativistic 7 factor of the * frame with respect to the lab frame, 
and similarly for the cm and projectile frames. (Notice that no particular particle is involved 
here. It’s just one frame with respect to another.) 

The quantity a* = ^ is defined as the speed of the * frame with respect to the lab divided 
by the speed of a particular particle in the * frame, and similarly for the cm and projectile frames. 

The quantity 07* = is defined as the speed of the * frame with respect to the lab di- 
vided by the speed of particle j in the * frame, and similarly for the cm and projectile frames. 

In general, xq refers to zeroth component of the 4-vector x where the components are written 
x 11 = (x°, x) = (.To, x). (The 0 is not a particle label. Particle labels are subscripts 1, 2, 3, 4, etc.) 

Mandelstam variable s = (pi -\- P 2 ) 2 ■ 
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Abstract 


In support of the development of 3-dimensional transport codes, this paper derives the 
relevant relativistic particle kinematic theory. Formulas are given for invariant, spec- 
tral and angular distributions in both the lab (spacecraft) and center of momentum 
frames, for collisions involving 2, 3 and n - body final states. 

1 Introduction 

One of the main new challenges in the development of radiation transport codes, such as 
HZETRN, is to provide for fully 3 - dimensional transport. This requires 3 - dimensional 
cross sections as input. Typically, these differential cross sections are calculated most easily in 
the center of momentum frame for nucleon - nucleon collisions, or the projectile frame for nu- 
cleus - nucleus reactions. However, transport codes require these cross sections in the lab frame. 
Even though there is extensive literature on transforming these cross sections from one frame 
to another, there is no reference that provides a complete discussion of this problem or includes 
special problems that need to be considered for radiation transport. One very common problem 
encountered in the literature is that lab differential cross sections are often written in terms of 
center of momentum or projectile variables. This is not useful for radiation transport because 
everything must be written in terms of lab variables. A major portion of the present paper is 
devoted to this problem. Significant complications arise because double - valued functions arise 
and these must be handled with extreme care. For example, two different angles in the center of 
momentum frame can correspond to a single angle in the lab frame. Another significant problem 
is that infinities can arise in the transformation of angular distributions to the lab frame. These 
infinities are rarely discussed in the literature. They are given significant attention in the present 
work. The aim of this paper is to provide a thorough analysis of relativistic transformations of 
differential cross sections, so that HZETRN can be upgraded for 3 - dimensional transport. 

This section discusses the formulation of differential cross sections in terms of scattering 
amplitudes, along with the various kinematic factors that are present. The types of differential 
cross sections that may be formed for 2 and 3 - body final states are also discussed. 

1.1 S-matrix and invariant amplitude 

The invariant amplitude is determined from the scattering matrix (often referred to as the S- 
rnatrix), and along with appropriate phase space factors, is used to determine the decay rate dr 
and cross section da differentials. Cross sections and decay rates are fundamentally calculated 
from a quantity called the scattering matrix or S-matrix, which is simply the time evolution 
operator for future and past times going to positive and negative infinity, S = U(+ oo,—oo). 
The Dyson expansion of the S-matrix is 

OO 

s = i + E 

71=1 

= Te~ i f d 


[ d 4 xid 4 x 2 ■ ■ ■ d 4 x n T[Hi(xi)Hi(x 2 ) ■ • • Hi(x n )] 

Til J 

ixH = dtH , (1) 
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where T is the time ordering operator, H is the Hamiltonian and 7i is the Hamiltonian density. 
Cross sections and decay rates are usually expressed in terms of an intermediate quantity called 
the invariant amplitude A4. which is related to the S-matrix. For the reaction 1 + 2 — ► 3 + 4, 
the relation is given by the Particle Data Group [1] and Griffiths [2] as 

(f\S\i) = {p3P4\S\pip 2 ) = 1 - iA/f(27r) 4 <5 4 (pi + p 2 - p 3 - Pi) , 1 , (2) 

y IL/2 & 

where p is the 4-momentum, E is the total energy, and the initial and final states are |i) = (pip 2 ) 
and |/) = (p 3 Pa)- With the above definition, the Feynman rules give —iM. In the above 
equation, the states are normalized according to 

{p'\p) = (2tt) 3 5 3 (p - p') . (3) 

Peskin defines the Transition matrix T (not to be confused with time ordering operator in the 
Dyson expansion) according to 


S = 1 + iT , 

and uses a different definition [3] (p. 104), 

(pi . . . Pn\iT\pAPB ) = iM{2-K) A 5 A {p A +PB - Zpf) , 


(4) 

(5) 


giving 

(f\S\i) = (pi . ..p n \S\pAPB) = 1 + iM{2ir) 4 5 A (p A +PB ~ £ Pf ) • (6) 

With the above Peskin definition, the Feynman rules give +iM. In the above equation, the 
states are normalized according to Peskin [3] (p. 23) 

(p|q) = 2£ , p(27r) 3 5 3 (p - q) . (7) 

For Dirac spinor states, Peskin [3] (p. 59) uses a similar normalization namely, 

(p, r|q, s) = 2 j F p (27t) 3 <5 3 ( P - p')^ . (8) 


Finally, note that Greiner uses a definition with different normalization for scalar and spinor 
states. See references [4] (p. 267), [5] (p. 221); 


(f\S\i) = iM{2TT) 4 5 4 (pi + p 2 

where the normalization factors are TV,; 
reference [4] (p. 267). 


-E^) n 


Ni 


n 


N- 


ti V 2^(2tt) 3 " y 2E' i (2ir) 3 ’ 

= 1 for scalar bosons and Ns = 2m. for fermions. 


(9) 

See 
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1.2 Phase space 

The phase space of a reaction contains all allowed momentum states of the final state particles. 
The phase space factor d4> n is used in expressions for dT and da. The decay rate T and total 
cross section a are obtained by integrating over phase space. 


1.2.1 Units and the invariant amplitude 

In this paper, natural units are used, with fi = c = 1, which means that both length and time 
have units of GeV” 1 and mass has units of GeV. Thus, cross sections a have units of GeV” 2 . 
The decay width T is related to the lifetime r of a particle by 
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(10) 


In natural units, the decay width has units of GeV giving the correct units of time, GeV” 1 , for 
the lifetime. The invariant amplitude A4 used below has units [2] (p. 200) of GeV 4 ” 11 where n is 
the number of external particles involved in a reaction. The number of external particles is equal 
to the number of external lines in the corresponding Feynman diagram. For example, when one 
particle decays into two particles, the number of external lines in the Feynman diagram is 3, 
and the units of M are GeV. For a 2 - body reaction producing two bodies in the final state, 
there are 4 external lines and so J\4 is dimensionless. 


1.2.2 dr and da in terms of the phase space factor 

The partial decay width (rate) of a particle of rest mass m decaying into n other particles is 
given by [1] 

dr = ^|M| 2 (2vr) 4 d4> n (P;p 1 , • • • ,p n ) , (11) 


or, in compact form as 

dr=^-|M| 2 (27T) 4 d$ n , ( 12 ) 

2m 

where M. is the invariant amplitude and S is a product [2] (p. 376) of statistical factors, 1 / j\ 
for each group of j identical particles in the final state, i.e. S = II a Jj. These statistical factors 
are discussed in references [3] (pp. 108, 151), [6] (p. 259), [7] (p. 110) , [8] (p. 200), [9] (p. 285) 
and [10] (p. 99). Gross [6] (p. 259) gives the most extensive discussion concerning the statistical 
factor versus limiting the regions of integration. The phase space factor is [1] 

n n dV 

d4> n (P;pi, •••,£>„) = £ 4 IP-X^)II / 9 , ( 13 ) 

i r=\ v /7 G mi 
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where n is the number of final state particles, and 


k 


P = Y,Pi > 

i= 1 


( 14 ) 


where k is the number of initial state particles. Equation (13) can be generated recursively [1], 
d$ n (P-,pi, ■ ■ ■ ,p n ) = d$j(q-,pi, ■ ■ ■ ,pj)d&n- j+ i(P-, q,p j+ 1 , • • • ,pn){2n) 3 dq 2 , (15) 


where 


3 

I 

1=1 


9 2 = (E^) 2 -IEp* 


i= 1 


Consider the particle reaction, 

1 + 2 — > 3 + 4 + • • • + (n + 2) , 

where the numbers label the particles. The microscopic differential cross section is [1] 

S 

da = —\M\ 2 (2ir) 4 d<S> n (p 1 +p 2 -,P3,---,Pn+2) 


S 
4 F 


\M\ 2 (2tt) 4 5 4 (pi +p-2 -P3-P4-..- Pn+2) 


d 3 ps d 3 p4 


d 3 Pn + 2 


(27t) 3 2E 3 (27t) 3 2E 4 (27r) 3 2E n+2 ’ 


or, in compact form [10] (p. 99), 


da = -^|Al| 2 (27r) 4 d4> n , 


(16) 


(17) 


(18) 


(19) 


where, again, S is a product [2] (p. 376) [10] (p. 99) of statistical factors, 1 / j\ for each group of 
j identical particles in the final state. The denominator flux factor F is defined as [10] (p. 98), 
[ 11 ] 


F = yf {pi.p 2 ) 2 ~ m{ml , 

where p\.p 2 denotes a 4-vector product. F can be written, 

+any frame — ^ ' J -^12 j 


where 


A jk = A (s,m 2 ,m 2 k ) 

= (s — m 2 — ml) 2 — Am 2 m\ , 


(20) 


(21) 


( 22 ) 


4 



with the Mandelstam variable. 


s = (pi +P 2 ) 2 • 


Note that 




fcj • 


Also, note that [10] (p. 98), 


Fi = m-2 |piz| , 


(23) 


(24) 


(25) 


and 


F c = Vs\pic\ = (E 1 + E 2 )\p lc \. (26) 

These equations agree with equation (5.113) in reference [11] and equations (4.52) and (4.53) in 
reference [10]. Note that F can be written in terms of the relative velocity v re i = \v\ — vz\. This 
method is used in references [3] (pp. 105, 106) and [12] (pp. 17, 18). In the above equation, the 
4-dinrensional delta function is 


<5 4 (P - Po) = S(E - E 0 )5 3 {p - p 0 ) , (27) 

or, for example, 

5 4 (pi + P2 ~ P3 ~ Pa) = S(Ei + E 2 - E 3 - E 4 )5 3 ( Pi + P2 - P3 - P4) , (28) 

where 5 3 (p — po) cancels a 3 - dimensional integral, with /( p) being an arbitrary function, as 


1 d 3 pf{ p) 5 3 { p - po) = /( Po) , 

(29) 

J d 3 p <5 3 (p - Po) = 1 • 

(30) 


1.3 Differential cross sections 

From equation (18), it can be seen that, in general, one will be able to form Lorentz invariant 
differential cross sections, such as d 3 a/(d 3 p/E), by bringing the relevant phase space factor to 
the left hand side. From such Lorentz invariant cross sections, one can form doubly differential 
cross sections such as, 


d 2 a , , d 3 a 
dEdn = |P| d 3 p/E ’ 


(31) 
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which follows from 


and 


and 


d 3 p 

dE 


p: 

E 


d\p\dfl , 


P‘ 


= E 2 - m 2 


2|p| d|p| = 2 E dE . 


(32) 


(33) 


(34) 


1.3.1 2 - body final state 

For a 2 - body final state, such doubly differential cross sections are meaningless, since E and 
6 are functions of each other. Therefore, either da/dE or da / dVt may be formed, but not 
d?a/dE dfl and therefore, not Ed 3 a/d 3 p. Consider the reaction, 

1 + 2 ->3 + 4, (35) 


where the numbers represent particles. From equation (18), the cross section is of the form, 

d 3 p 3 d 3 p 4 


d <7 rx\M\ 5{Ei + E 2 -E 3 - Ei)S ( Pl + P2 - P3 - Pi) (2t)32£ . 3 (2n)32Ea ■ 

The delta function, < 5 3 (pi + P2 — P3 — P4), cancels the integral over d 3 p 3 , resulting in, 

1 d 3 p 4 


(36) 


da oc \M\ 2 5{E l + E 2 -E 3 - Ea) 


(: 2it) 3 2E 3 {2tt) 3 2E 4 


(37) 


However, since d 3 p 4 = p\dp 4 dQ . 4 oc dE 4 d^l 4 , and the energy delta function cancels the integral 
over dE 4 , the result is 


da oc \M\ 2 dQ 4 . (38) 

Thus, for a 2 - body final state, one can form the single differential cross section da/dVL 4 , or 
alternatively da/dE 4 , because energy is a function of angle. 

Finally, note that for a 2 - body final state, the angular distribution da/dVL can be formed 
in both the lab and cm frames. However, the spectral distributions da/dE and da / dT , where 
T is kinetic energy, can only be formed in the lab frame and not in the cm frame. The reason is 
that both T and E = m + T are fixed in the cm frame, but vary with angle in the lab frame. An 
exception to the above statement is if the produced particle has a variable mass. An example is 
the A particle. (See reference [13].) However, a variable mass implies that the particle decays, 
which means that ultimately the final state is not a 2 - body state. 
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1.3.2 3 - body final state 

Now, consider the 3 - body final state, 

1 + 2— *3 + 4 + 5. (39) 



1.3.3 Unit conversion to mb 

As shown above, the units of da/ctt are GeV~ 4 . These units are converted to nrb/GeV 2 by the 
conversion factor GeV~ 2 = 0.389379 mb. 


1.4 Proof that the differential cross section Ed 3 cr/d 3 p is Lorentz invariant 

A proof of the Lorentz invariance of Ed 3 cr/d 3 p can also be found in reference [14] (p. 52). 
Differentiating p 2 = E 2 — p 2 gives, 

d 4 p = dE d 3 p . (45) 

Expanding 5(p 2 — m 2 )d 4 p gives, 

5{p 2 — m 2 )d 4 p = 5(E 2 — p 2 — m 2 )dE d 3 p . (46) 
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Applying the formula 


S(f(E)) = Y,j^- ) S(E-E i ), (47) 

with f(E) = (E + a/p 2 + m 2 )(E — v^P 2 " + m 2 ), and keeping only the positive E root, gives 

6(p 2 — m 2 )d 4: p = S(E — \fp ^\-m?)dE—^ . (48) 

v 2 E 

Integrating over dE cancels the delta function, leaving 

2^ = <^'(P 2 - • (49) 

Since p 2 and m 2 are Lorentz invariant, the right hand side is invariant, and therefore d 3 p/2E 
is invariant. The total cross section, a, is invariant since da is transverse to the direction of 
motion z, and transforms like an area. Consequently, Ed 3 a/d 3 p is invariant. 

As an alternate proof, expand the momentum integration measure, 

d 3 p = dp z dp x dpy . (50) 

Here, dp x and dp y are Lorentz invariant, being transverse to the direction of motion z. Since a 
is Lorentz invariant, it remains to show that Jd is invariant. Rapidity y, is defined through 

E = hit coshy and p z = mx sinh y , (51) 


and 


rriT 


2 i 2,2 

m + px + p y 
E 2 - pz 2 , 


(52) 


which gives 


dp z 

dy 


rriT cosh y = E . 


(53) 


Thus, 


E _ 1 
dp z dy ' 


(54) 


Under Lorentz transformations, rapidities add, so that y = y' + yup , where yaF is the rapidity 
of a certain reference frame, y' is the rapidity in that reference frame, and y is the transformed 
rapidity. Thus, dy = dy ' , showing that dy = dp z /E is invariant under Lorentz transformations. 
Thus, Ed 3 a/d 3 p is Lorentz invariant. 



1.5 Denominator flux factor in any frame 

The denominator flux factor F is defined in reference [10] (p. 98) as 


This can be written 


where [14] (p. 23), [15] 


F = \J (pi+ 2) 2 ~ rn\m% . 


Fa.ny frame \/Al2 j 


A (x,y,z) = (x - y - z) 2 - 4yz 

= x 2 - 2(y + z)x + (y - z) 2 


= x 2 + y 2 + z 2 — 2 yx — 2 zx — 2 yz , 


and 


A ij = A (s, m 2 ,m 2 ) = (s — m 2 - m 2 ) 2 - 4 m 2 m 2 . 

Now calculate s as, 

s = {pi + P2)' 2 = pi + pi + 2 pi ■ P2 

= m\ + to 2 + 2 {E\E 2 — pi • P 2 ) in any frame . 


Then, 

A 12 = A(s, ml, m^) = (s — rn\ — m\) 2 — 4 m\m\ 

= 4[(E\E-2 - pi • p 2 ) 2 - to^toI] 
= 4[(pi • p 2 ) 2 - rn\ml\ . 


( 55 ) 


(56) 


(57) 


(58) 


(59) 


(60) 


1.6 Denominator flux factor in the lab frame 

For the reaction 1 + 2— ► 3 + 4 + • • •, let the rest frame of particle 2 be the lab (target) frame. 
Then, P 2 = 0, and 

(p\ ■ P2) 2 ~ m\m\ = [E x E 2 - pi • p 2 ) 2 - to 2 to| 

= E\E 2 2 — to 2 to 2 (because P 2 = 0) 

= (Piioft + ™\)m 2 2 - m\m 2 2 

= "I 2 P llab , (61) 


giving 


Flab = rn2\Plllab 


(62) 


9 



2 Decay into two particles 

This section is concerned with the decay of one particle into two particles. The fundamental 
quantity is the decay width. Formulas are given that enable this to be calculated from the 
invariant amplitude. The kinematics of the final state particles is also discussed. Finally, the 
decay width is related to the cross section for the reverse process in which two particles react to 
form a single particle. 

2.1 Decay width 

Consider the reaction, 


Ni + N 2 - N! + A - IVi + vr + iv 3 , (63) 

where N t are nucleons (either neutron or proton), A is the A particle and n denotes a pion. 
Pilkuhn [16] (pp. 37, 38) and de Wit [7] (p. 113) show that the cross section for this two step 
process can be approximately written as, 

rf A — » 7 r + N‘*) 

a(N 3 + N 2 -> JVi + vr + 1V 3 ) « a(Ni + N 2 -y Nx + A)— T • (64) 

In the above, F is the total decay width for the process A — ► anything. We seek expressions for 
2 - body cross sections, as well as the decay of a single body into 2 bodies. Consider the decay, 


3 — > 1 + 2 , 

where the numbers denote particles. The 2-particle decay width is given by 

dr = ^|M| 2 (27r) 4 d$ 2 (P; Pl ,p 2 ) , 

2ru. 


with m = m3. The phase space factor is 


d$ 2 (P',Pi,P2) = S 4 (p 3 -pi-p 2 ) 


d 3 pi d 3 p2 
{2tt) 3 2E 1 {2tt) 3 2E 2 


— 5(E 3 — E\ — Ei 


,/ \ d 3 pi 


= 5(E 3 — E\ — Ei 


eie 2 

|pi| 2 d|pi|dPi 


eie' 2 


(65) 


( 66 ) 


(67) 


with 


e = {2tt) 3 2E , (68) 

where it is understood that 

P 2 = P3 Pi ) (69) 
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so that in the above equation, with |pi| 2 = p 2 , 


Ei = y/p\ + m\ (70) 

E' 2 = p 3 - pi) 2 +mj . (71) 

Two equivalent methods to eliminate the 5 function are now given. These are the Ho-Kirn phase 
space method [10] and the Griffiths phase space method [2]. 

2.1.1 Ho-Kim phase space method 

Ho-Kirn [10] eliminates the 8 function by making a change of variables as follows. Simplify 
equation (67) by writing 


+ (72) 

so that the integral over d(E i + E 2 ) will cancel the deltal function 8 ( E 3 — E\ — E 2 ). Ultimately, 
the quantity Efpif 3 ' 1 will be evaluated. The quantity, d(E\ + E! 2 ), is used rather than dE \ , 
because both E\ and E 2 contain the term |pi|. The phase space factor becomes, (with the energy 
5 function eliminated), 


d<r 2 (P;pi,p 2 ) 


|Pi| 2 ^|Pi| 

eie' 2 d(Ei + E 2 ) 


dVL\ 


(in general), 


giving the partial width as, 


dF 


S 

327r 2 ?n 


\M\ 2 


|pi| 2 ^|Pi| 

EiE' 2 d(Ei + E 2 ) 


drii 


(in general). 


(73) 


(74) 


2.1.2 Griffiths phase space method 

Griffiths [2] eliminates the 8 function by making use of the formula, 


%(x)] 


n 


E 


d{x - Xj) 

W(Xj) | 


(75) 


where the sum over j represents the sum over the zeroes xj of the function g{x). The term 
g'(xi) is the derivative of the function g(x), evaluated at the zero, Xj. In the delta function, 
8(E 3 - Ei - E' 2 ) = <5[ff(x)], we set, 

g{x)=g{ |pi|) = E 3 -Ei-E 2 

= E 3 - y/pl + m\ - \J(P3-Pi) 2 + m 2 1 . (76) 
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The zeroes, pij, are obtained by setting 


s(|Pi|) = E 3 - ypf+mf - sj (p 3 - pi ) 2 + rrij = 0. 


(77) 


The derivative is 


_ rfg(|pi|) _ 4^1 + ^ 2 ) 

d|pi| d|pi| 


(78) 


so that 


S(E 3 -E x - ^ 2 )d| Pl 


5(|pi| - |pii|)d|pi 
d(Ei+E |) 

«2|pi| 


<^|Pl| 

d(E 1 + E! 2 ) ’ 


(79) 


which is the same result obtained with the Ho-Kirn phase space method. The phase space factor 
and the decay width are evaluated in a particular frame by evaluating y in that frame. 

The cm frame, in which the decaying particle is at rest, is the most logical choice of frame. 
Thus, 


P3 = 0 

and E! 2 = ^\pi\ 2 + m%, giving, 


and 


d(E x + E' 2 ) 

d|Pi| 


= Pi 


Pi = -P 2 


, E! + Ek 


E\E’ 


Since E\ + E 2 = £3 = m and |pi| = |p 2 1 = |p/| in the cm frame 

12 d| Pl ' 


d$ 2 {P',Pi,P2) = 


Pi I 


eie^ d{E\ + E 2 ) 
I P 1 1 2 E\ E 2 


Trdtt 1 


eie 2 m 
|Pi| 

4(27t) 6 m 


‘-dQ 1 


d£l± . 


(80) 


(81) 


(82) 


Assuming no angular dependence in At, the integral / dTl 1 = 47 T, gives the width in the cm 
frame, 


U = ^l^l 2 <2») 4 «2(P;Pi,P2) = ■ (83) 

This is only valid if M has no angular dependence. This agrees with references [2] (equations 
6.30, 6.32) and [10] (p. 103) (with dQ = 4 tt). The units can now be verified. . There are 3 
external lines in the Feynman diagram, so that M has units of GeV. The units of M 2 cancel m 2 
in the denominator, leaving T with the correct units of GeV. Using the result from the following 
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section. 


| P 2 1 = |p/| = — — - ymf + mf + m\ — 2mlml — 2m\rn^ — 2mlml , 

which is equivalent to solving equation (77) for the zeroes. The decay rate is written 


T c = ^ 3 \M \ 2 y 777,3 + + 777,3 — 2 mlml — 2 rn 2 ml — 2m2?n 3 . 


(84) 

(85) 


with m = 777 , 3 . Again, this is only valid if M has no angular dependence. This agrees with 
reference [10] (p. 103). 


2.2 Decay energies and momenta 

The energy E 2 , and 3- momentum P 2 , of the decay, 

3 — > 1 + 2 (86) 

may be expressed in terms of the masses of the three particles. Conservation of 4-momentum 
gives 

p\ = (P 1 + P 2) 2 

= m l = P 1 +P 2 + 2 Pi-P 2 = ml + m\ + 2{E x E 2 - P 1 .P 2 ) 

= ml + m| + 2{E\E 2 + p 2 ) (because pi = — P 2 in the cm frame) 

= ml + ml + 2{E\E 2 + E% — m|) 

= T m-j + 2E 2 {E\ + E'2 ) — 2 ml = ml — ml 2E 2 {E\ + E 2 ) . (87) 


Also, 


Pi = {P1+P2) 2 

= ml = {E x + E 2 ) 2 - (pi + p 2 ) 2 , 


( 88 ) 


giving, 


7773 = Ei + E 2 , 


(89) 


since Pi + P 2 = 0 in the cm frame. Substituting equation (89) into equation (87) leads to 


=$■ e 2 


ml - ml + 2E 2 ms , 
ml + rn 3 ~ m l 
2ms 


(90) 

(91) 


The magnitude of the 3-momentum is given by 
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2.3 Two - body reaction producing one body 

We can relate the decay width of the process, 

3 — > 1 + 2 , 


to the cross section for 


(92) 

(93) 


1 + 2 -> 3 . 

The cross section is given by equation (19), 

da = -^\M\ 2 (2ir) 4 d<5> n (pi+p 2 -,P3) 

= ^\M\ 2 {2ir) 4 5 4 {pi +P 2 - p 3 ) (many frame), 

= j r ^j-^|Al| 2 (5(-E , i + E 2 — -E 3 )i5 3 (pi + P 2 — P 3 )-yy^ (cm frame), (94) 

4|pi|V s e 3 

using equation (26). In the cm frame, we have E 3 = m 3 = m = y/s and pi = — P 2 , finally giving 
the cross section, (with 5 = 1 for the above reaction 1 + 2 — > 3), 

da = I \M\ 2 8{Ei + E 2 - E 3 ) , (95) 

4m/ 1 pi | 

with E\ = \Jy>\ + m\ and E 2 = \Jp>\ + rn\. We now relate this cross section to the decay width. 
Substituting |Af| 2 from equation (85), into equation (95), gives the relation between decay width 
and cross section as, 

27T 2 

da = — To r S(Ei + E 2 - E 3 ) . (96) 

5|pip 

The units of 5{E\ + E 2 — E 3 ) are GeV -1 , cancelling the GeV units of T, so that the units of a 
come out correctly as GeV -2 . 

Sometimes the delta function is written a little differently, in terms of 5(s — m 2 ). In the cm 
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( 97 ) 


frame, y/s = E\ + E 2 = W 2 . Thus S(s — m 2 ) = 5{W 2 — m 2 ). Using the result [2], 

5(x 2 - a 2 ) = — ^j[(S(x — o) + 5(x + a)] , 


gives 

S(s -m 2 ) = 5[{Ei + E 2 ) 2 - m 2 } = ^ - [5(Ei + E 2 - m) + 6 {E 1 + E 2 + m)} . (98) 

Zi II h 

However, we cannot have E\ + E 2 = —m, so that the second term does not contribute, [16] (p. 
15), [3] (p. 151), giving only 


5(s - m 2 ) = -^<5(£i + E 2 - m ) . 

2m 

Thus, the relation between the decay width and the cross section can also be written 

47r 2 m 


da = 


«S|Pi| 


T S(s — m z ) . 


(99) 


( 100 ) 


Suppose the decay products are both massless (e.g. H — > 77 ). Recall that in the cm frame, 
the momentum of the decay products is |pi| = |p 2 1 - Also, the mass of the decaying particle is 
E 3 = m 3 = m. By energy conservation (with E = |p| for photons), we have 

m = 2 |p | . ( 101 ) 

Substituting |p| = m/2 gives, 

da = — — T 5(s - m 2 ) . (102) 

Sm 

Now, if the particles have spin, then it must be averaged over the initial spin states S\ and S 2 
and summed over the final spin state S3 [7] (p. 114). This brings in an extra factor 


da 


2S 3 + 1 16tt 2 

(25i + 1)(25 2 + 1) Sm [S 



(103) 


which is our final result. This gives the results of several other authors. For example, Norbury 
[17] considered the decay of a Higgs boson into two photons, 


H -> 7 + 7 , 


where the statistical factor is S = 1 / 2 !, because there are 2 particles in the final state. The spin 
of the Higgs is Sh = S 3 = 0, and the spin of the photons are Si = S 2 = S 7 = 1. However, 
the factor, 2g 1 +1 = 1/3, is only correct for massive spin 1 states. Photons loose one degree of 


15 



freedom, meaning that the factor should actually be 1/2. Substituting gives 

o 2 

da = T S(s — m 2 ) , (104) 

m 

in agreement with [17]. Another example is the decay of a vector boson B into an electron- 
positron pair [3] (p. 151), 


B -> e + + e~ . 


Here, 5 = 1 because there are no identical particles in the final state. Also, the spins are 
Sa = Sb = 1, and Si = S e + = 1/2, and S 2 = S e - = 1/2. Substituting, gives 


da=- 2n2 ^ x ^ 


m 


T 5(s — m z ) , 


(105) 


in agreement with [3] (p. 151). Finally, note that Vidovic et al [18] have an extra factor of 2 in 
their equation (46). However, this results from the definition of their equation (21). 


2.3.1 Calculation using de Wit and Smith formula 

De Wit and Smith [7] (p. 114) have also derived a relation between the cross section and decay 
width that is more general than our result. In this section, we show that the de Wit and Smith 
formula, reduces exactly to our result. They consider the reaction 


a A b — ■> A — > A , 

where A has decayed into X, i.e. A — > X. Their result is 

T{A^ X)T{A->a + b) 


a(a + b — > X) = 


I67 Tm\ 


(s\ — m\) 2 + m^T' 2 


where S is the same statistical factor as the present work and X(m\, m 2 , ?n 2 ) 
as used herein, namely 


(106) 

(107) 

is also the same 


\{m\,m 2 a ,ml) = 4|pi c | 2 s = m\ , 


(108) 


where we have used equation (101). The above result is valid for arbitrary width. Now, examine 
the narrow width approximation. A narrow width corresponds to a large lifetime. In other words, 
the narrow width approximation has as its limit the decaying particle not decaying. That is, the 
decay A — > X does not proceed. Thus, the symbol X in the above equations should be replaced 
with A. In the narrow width approximation, one has [7] (p. 112) 


1 m A T{A -»• X) 

7 r (s\ — m\) 2 + m\T' 2 


5(s a - m\) . 


(109) 
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Substituting equations (108) and (109) into equation (107) gives 


a (a + b — > A) = — r(A — > a + b) 5 {sa — m 2 A ) , (110) 

orriA 

which is the same as our result in equation (103). 

3 Non - relativistic two - body kinematics 

Some key results in this paper will be in finding the relation between the lab angle and the lab 
momenta for 2 - body final state particles. An important outcome is that for a given lab angle, 
there will be two corresponding momentum values, which result from two angles in the center 
of momentum or projectile frames. See equations (146) and (147). This is a crucial point for 3 - 
dimensional transport codes, because for a given angle in the lab (spacecraft) frame, there will be 
two particles arriving from the center of momentum or projectile frames with different energies. 
The differential cross section for the production of each particle must be added incoherently to 
form the differential cross section in the lab frame. This section examines these results from 
a non - relativistic point of view. This will enable easier understanding of the corresponding 
relativistic results, which will be derived later. 

3.1 Reactions between elementary particles 

First, examine collisions from a more general point of view. Consider the collision of elementary 
particles given in the reaction 


1 + 2— >3 + 4, (111) 

where the numbers represent the particles involved in the collision. At energies comparable to the 
masses of the particles, relativistic kinematics is needed to accurately analyze the collisions. This 
section considers low energy collisions, and will use non-relativistic kinematics. By conservation 
of mass 


m\ + m 2 = m -3 + m 4 . 


( 112 ) 


Momentum is always conserved, 


Pl + P2 = P3 + P4 (113) 

We don’t need initial i and final / subscripts now because particles 1 and 2 are automatically in 
the initial state and particles 3 and 4 are automatically in the final state. For elastic collisions, 
kinetic energy will be conserved, 


Ti + t 2 = r 3 + t 4 . 


(114) 
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Remember that the kinetic energy and momentum can be expressed in terms of each other via 


T = 



(115) 


which is the non-relativistic version of kinetic energy. Also, note that we will often use the 
identity 


P 


2 



(116) 


where |p| is the magnitude of the momentum vector p. In particle reactions like (111), particles 1 
and 2 are usually produced in a particle accelerator and are made to collide with each other. The 
mass and energy or momentum of particles 1 and 2 is always known. Thus, mi, m 2 , Tj, T 2 , Pi, P 2 
are always known. 

In a fixed target accelerator , the target particle 2 is at rest and the projectile particle 1 is 
fired into it. The collection of projectile particles is often called the beam. The collection of 
target particles is often just called the target. Thus, the moving accelerator beam is fired into a 
target at rest. The whole accelerator complex is called a laboratory or simply a lab. Thus, the 
target particle is not moving in the lab. We often call this reference frame in which particle 2 is 
at rest simply the lab frame or the target frame. So, in a fixed target accelerator we have 


P 2 = 0 and T 2 = 0 . 


(117) 


In a colliding beam accelerator , often called a collider , both particles 1 and 2 are moving in 
opposite directions and they are fired into each other. However, again we will always know the 
masses and energies of both of these particles. 

In either case, particles 3 and 4 are the reaction products and their properties are studied 
in great detail. These properties include their mass, charge, spin, and other quantum numbers, 
and also their kinematic properties such as their energy or momentum and also the angle from 
which they emerge from the collision. 


3.1.1 Mandelstam variables 


If we write down the x and y components of the momentum conservation equation (113), we can 
then deduce the angles with which particles 3 and 4 emerge. However, this can involve quite a 
lot of algebra. An easier method is to use the so-called Mandelstam variables, which are defined 
in their non - relativistic form as 


(Pi + P2) 2 = (P3 + P4) 2 , 

(118) 

(Pi - P3) 2 = (P4 - P2) 2 , 

(119) 

(Pi - P4) 2 = (P3 - P2) 2 • 

(120) 


These are just the usual relativistic Mandelstam variables with 4-vectors replaced by 3-vectors. 
Note that s is related to the total energy of the system and t and u are momentum transfers. 


18 



The square of the sum of any two 3- vectors is 


(A + B) 2 


A 2 + B 2 + 2A • B , 

A 2 + B 2 + 2| AllBl cos 0 


( 121 ) 


where 0 is the angle between the vectors, as seen in figure 1. The Mandelstam variables can 
then can be written 


s 

t 

u 


Pi + p| + 2 |Pl||P2| COS 012 
P 1 +P 3 - 2|pi||p 3 | COS 013 
Pi +P4 - 2 1 Pi I |p 4 1 COS 014 


P3 T P4 T 2 1 p 3 1 1 p 4 1 cos 0 3 4 , 
P2 + P4 - 2|p 2 | |p4| COS 024 , 
P 2 + P 3 -2|P2||P3|COS023 , 


(122) 

(123) 

(124) 


where 9jk is the angle between particle j and k, as shown in figure 2. 



Figure 1: Scalar product between two vectors A and B is given by A.B = AB cos 6, where A 
and B are the vector magnitudes. Note that the angle 0 is defined as the angle between the 
vectors arranged tail to tail. 



Figure 2: The angle 9jk is the angle between the vectors pj and pk connected tail to tail as in 
figure 1. The scalar product is given by Pj-Pfc = | P j 1 1 P fc | cos Qjk- 


3.1.2 Lab frame 

As discussed above, the lab frame is the frame in which the target particle 2 is at rest as in 
equation (117). With p 2 ; = T 2 / = 0, our conservation equations become 

Pk = P3 + P4 (125) 
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and 


and if the collision is elastic, 


mi + m2 = m3 + mi , 


Tu = T 3 + T 4 . 


The momentum conservation equation can be written in terms of its components as 


(126) 


(127) 


x component : 

mivu = 777 , 31/3 cos 0 i 3 + 7 / 741/4 cos 0 i 4 , 

(128) 

y component : 

0 = 77731/3 sin 0 i 3 — 7 / 741/4 sin 0 i 4 , 

(129) 


or as 


x component : 

Pu| = P3 COS 013 + P4 COS 014 , 

(130) 

y component : 

0 = P 3 sin 0 i 3 — p 4 sin 0 i 4 . 

(131) 


Given that the target particle 2 is at rest, then it is meaningless to define an angle relative to 
particle 2, since angles are defined between lines and not a line and a point. Thus, in the frame 
in which particle 2 is at rest, there will be no angle appearing that contains particle 2. That is, 
there will be no angle 0j2- Setting p 2 / = 0 gives the lab frame Mandelstam variables 

S = Pi/ = p! + P4 + 2 Ip 3 ||P4|cos034 , (132) 

t = Pi + p| - 2|P! I |p 3 | cos 6> 13 = ph , (133) 

u = P? + P 4 - 2|pi||p4|cos0 14 = p| • (134) 

The portions of these equations that do not have a l subscript are true in any frame. The angles 
are shown in figures 1-4. Note that 


034 = 013 + 014 • 


(135) 


3.1.3 Proof of the relation t = (pi — P 3) 2 = |P 4/| 2 

Using the x and y components of the momentum conservation equation, one can prove the 
Mandelstam variable relation 

t = (pi — P 3) 2 = |piz | 2 , in the lab frame where p 2 / = 0. (136) 

This statement will now be proved. Conservation of momentum is 

P 1 + P 2 = P3 + P4- (137) 
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1 



Figure 3 : Definition of angles in the lab frame. 0 m/ is the recoil angle of the target. The angle 
between the final state particles is 6*34 = #13 + 0 m- Note that the angle 6 is defined as the angle 
between the vectors arranged tail to tail. 



Figure 4 : Definition of angles in the cm frame. Note that 0i3 C = 024 c and 0 Mc = 023 c- The angle 
between the final state particles is 034 = 0i3 + 0i4. In the cm frame, 034 = 7 r. Note that the 
angle 0 is defined as the angle between the vectors arranged tail to tail. 


All quantities are now to be considered as lab frame quantities. The components are 


x component : |pi| = | P3 1 cos 0i3 + | P4 1 cos 0i4 , 

y component : 0 = | P3 1 sin 0i3 — | P4 1 sin 0i4 . 


Now eliminate | p4 1 and 0 m. The y component equation gives 


sin 0i4 


P 3 | sin 0 j 3 

|P4 1 


Use cos 0 


±71 


sin 2 0 and substitute into the x component, 


Pi 


p 3 | cos 013 ± |p 4 



I 12 • 2 

|P3T sm 

IpII 


013 


( 138 ) 

( 139 ) 


(140) 
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(141) 


p 3 | COS $13 ± 



2 sin 2 0 i 3 . 


This gives 


(|Pi I — |P 3 | cos 6»i 3 ) 2 = |p 4 | 2 - |p 3 | 2 sin 2 0 i 3 

= |pi | 2 - 2 |pi||p 3 |cos 0 i 3 + |p 3 | 2 cos 2 0i 3 , (142) 

so that 

I P 4 1 2 = |Pi | 2 - 2|pi||p 3 |cos0i 3 + |p 3 | 2 = (pi - p 3 ) 2 . (143) 

Note that when solving collision problems with the x and y components of momentum conserva- 
tion, one always ends up with the Mandelstam variable relation anyway, if angles are of interest. 
Thus, it is good practice to simply start with Mandelstam variables when angles are involved. 
Squaring the quantity, 


(Pi - P 3) 2 = t , (144) 

will directly give the angle $ 13 . Similarly, squaring the quantity, 

(Pi - P 4) 2 = u , (145) 

will directly give the angle $ 14 . Now, a very interesting result is that in the lab frame, for 
reaction (111), the speed or energy of particle 3 is given in terms of the angle at which particle 3 
emerges, namely 0 43 . Similarly, the speed or energy of particle 4 is given in terms of the angle at 
which particle 4 emerges, namely 0 44 . In other words, the angles of the final particles determines 
their energy and vice-versa. 

3.2 Relations between momentum and angle of scattered particle 
3.2.1 Center of momentum frame 

Consider the reaction, 


1 + 2— >3 + 4 

in the cm frame. As will be shown, the momentum of particle 3 or 4 is not determined by the 
angle at which it emerges from the collision. 

3.2.2 Lab frame 

In the lab frame, the momentum of particle 3 or 4 does depend on the angle at which it emerges 
from the collision. This is shown below. Again, consider the elastic reaction 

1 + 2 — > 3 + 4 , 
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where the target particle 2 is at rest in the lab and the projectile particle 1 is fired into the 
stationary target. Assume that we know the value of the momentum of particle 1 to be pi/. 
One can prove the following useful results. 


(a) The momentum of particles 3 and 4 can be expressed in terms of their scattering angles 
$13/ and $14/ , where dm is the angle between the momentum of particle i and the incident 
beam direction of particle 1 , as measured in the lab frame. One can show that 


mil = Pn 


m 3 (cos $i 3 z ± yiflf - sin 2 $^ ) 


777-3 + 7774 


( 146 ) 


and 


|P 4 Z| = |P1Z| 


7774 (cos $ 14 / ± - sin 2 $u/ ) 


( 147 ) 


7773 + 777-4 

Note that |p4/| is the same as |p3/|, except for the interchange of the particle labels 3 <-> 4 . 


(b) We will evaluate |p3/| and |p4/| for an elastic billiard ball collision when mi = m3 & m2 = 
7774 and also when mi = 7774 & 7772 = m3. 


(c) Notice that the solutions for |p3/| and |p4/| are double valued. For the case of an elastic 
billiard ball collision when mi = 777.3 & m 2 = 7774, we will interpret the two solutions. 

(d) We shall evaluate |p3/($i3/ = 0,7r, |)| and |p4/($i4/ = 0,7r, §)|. We shall also work out 
these expressions when mi = m3 & 7772 = 7774. 

(e) We shall interpret the solutions |p3/($i3/ = 7r/2) | and IP47 ($147 = 7r/2) | when mi = 
m3 & 777-2 = 7774. 

(f) Consider a 1 -dimensional collision where $13/ = 0 or 7r and $14/ = 0 or 7r. We have already 

evaluated IP37 ($137 = 0 , 7r) | , and |p4/($i4/ = 0,7r)| for an elastic billiard ball collision when 
777-1 = 7773 & m2 = 7774. We will interpret the results for the three special cases i) 

mi >> 777-2, ii) 7771 = 7772 , and hi) 777 l << m-2- 


We shall now prove the statements given above. All quantities are assumed to be lab frame 
quantities. 


(a) Conservation of momentum, mass and energy is 

( 148 ) 

( 149 ) 

( 150 ) 

We use the Mandelstam variable, 

t = (pi -p 3 ) 2 = (P 4 -P 2) 2 


Pi + P2 = P 3 + P 4 , 
?77l + m2 = m3 + 7774 , 

T l + T 2 = t 3 + t 4 . 
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(151) 


= |pi | 2 + |P3| 2 - 2|pi I |p 3 | COS 6>13 


= P4 


in the lab frame where P 2 = 0 . 


Note that 


i 1 2 _ 2 

P =P 


Conservation of kinetic energy can be written (with T 2 = 0), 


ipir 

2mi 


I P3 1 2 + | P4 1 2 
2 m .3 27774 ’ 


which gives 


2 ( IpiI 2 IP 3I 2 ' 

P 4 = 7774 

\ m 1 7773 


so that the Mandelstam variable t becomes 


1 1 2 1 1 1 2 01 11 1 n 771-4 2 7714 2 

Pi + P3 - 2 Pl P3 cos 013 = P! P 3 

mi m 3 


from which we get the quadratic equation for | P 3 1 , 


m .4 


P 3 1 + — - 2|pi| |p 3 | cos 6 »i 3 + |pi|- 1 = 0 . 


777.3 

A general quadratic equation can be written 

ax 2 + bx + c = 0 , 

with solution 

—b ± Vb 2 — 4ac 


7774 \ 


7771 / 


X = 


2 a 


For our equation, we have 


b 2 — 4 ac = 


4| P i| 2 cos 2 0 i 3 - 4 1 p 1 1 2 [l - (l + 

V mi) \ m 3 ) 


4|Pi| 2 

4 Ipi | 2 

4|Pi | 2 


cos^^^-fl-^ 4 ) { 1 + " 14 

7771 / V 7773 / J 


2 a , , 7714 7714 , 777| 

COS 013 - 1 H 1 

mi m 3 m±m 3 

— sin 2 013 H (77737774 — 77717774 + m 4 ) 

mim3 


(152) 

(153) 

(154) 

(155) 

(156) 

(157) 

(158) 
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4 1 Pi 1 2 

7774 

.mi 7773 

4|pi| 2 1 

/ 77727774 

\777im.3 


(m3 — mi + m.4) — sin 2 $13 
■ - sin 2 0i3 


( 159 ) 


where we have used conservation of mass, m3 — mi + mi = m2, in the last equation. Thus, 
the solution to the quadratic equation is 


P3 = 


2 1 pi | cos 013 ± 2 |pi|y / '(^| - sin 2 0 i3 


2 (! + ^) 


= Pi 


m 3 (cos 0i3 ± yiflf - sin 2 0^ ) 


777-3 + rn 4 


( 160 ) 


which is the desired result. We could work everything out again to determine |p4 1 . How- 
ever, rather than going through all the algebra again, it should be obvious that our answer 
is the same as that obtained previously except that subscript 3 is replaced by 4 . Thus, we 
immediately deduce that 


P4 = iPi 


7774 (cos 0i4 ± -sin 2 0^ ) 


( 161 ) 


m3 + m 4 

The Mandelstam variable that will appear is 

u = (pi — P4) 2 = | P3 1 2 , in the lab frame where P2 = 0. ( 162 ) 

(b) If 777 i = m 3 & 7772 = m 4 (elastic billiard ball collision), the above results reduce to 


P3 = Pi 


P4 = Pi 


= Pi 


7771 (^COS 013 ± y ^ - sim 013 

mi + m2 

m 2 (COS 014 ± COS 014) 

7771 + m2 

2 m -2 cos 0i4 


mi + m2 


or 


0 . 


( 163 ) 


If 7?7i = 7774 & m2 = m3, we get 


P3 = IPil 


= Pi 


m2 (cos 013 ± COS 013 
777i + m 2 

2m2 cos 0i3 
777i + m 2 


or 0, 


( 164 ) 
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and 


P4 


Pi 


m i 


^COS 6*14 ± 



mi + m2 


(c) For mi = m 3 & m 2 = 111 . 4 , we had 


P3 


Pi 


mi 


^COS 6*13 ± 



mi + m2 


and 


2 m 9 cos 0 i 4 

P 4 = Pi ; or 0. 

m 1 + m 2 


(165) 


(166) 


(167) 


The identification, mi = m 3 & m 2 = 111 , 4 , means that particle 3 is the scattered projectile 
and particle 4 is the recoil target. First, we interpret the | P 4 I = 0 solution. This implies 
that P 4 = 0. We already have P 2 = 0 in the lab frame. Thus, conservation of momentum 
Pi + P 2 = P 3 + P 4 implies that 


Pi = P3 , 


(168) 


and together with 


P2 = P4 = 0 , 


(169) 


this means that the projectile misses the target. Thus, the target remains undisturbed 
(p 4 = 0 ) and the projectile continues on with its original momentum (p 3 = pi) as if 
nothing happened. Obviously, the other solution | P 4 1 = |pi| 2r ^ c °^ 4 signifies a real 
collision, not a miss. Now we interpret the |p 3 | solutions. Firstly, we consider whether 
there is a miss solution as above. The miss solution corresponded to pi = P 3 which implies 
that 0 i 3 = 0. Thus, if $13 7 ^ 0, a real collision has occurred and the momentum |p 3 1 is 
simply a double valued function of the scattering angle $ 13 , and no further interpretation 
is warranted. We look at the $13 = 0 solution. This gives 


P3| = IP 1 


(170) 


or 


. . . mi — m2 

P3 = Pi : • 

mi + m2 


(171) 


Clearly |p 3 1 = |pi| corresponds to the miss solution and | P 3 1 = |pi| corresponds to 

the real collision. 
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(d) Using equations (146) and (147), we obtain 


|P3(0i3 = O)| = |pi| 

= IPil 

|P3(^13 = 7t)| = “IPil 

= “IPil 

|P3(6 > 13 = vr/ 2 )| = ±|pi| 

= ±|Pi| 

|p 4 (6»i4 = 0)| = |pi| 

= IPil 

| P4 (6»14 = 7r)| = -|pi| 

= -IPil 


m 3 + m .4 

mi ± m 2 
mi + m 2 ’ 

m 3 

m 3 + m. 4 

mi =F m-2 


1 ± 


mim 3/ 

when mi = m 3 & m 2 = m .4 , 


/mamA 
V mim 3 / 

when mi = m 3 & m 2 = m 4 , 


mi + m. 2 

m 3 /m 2 m 4 “ 

m 3 + m .4 V mim 3 

m 2 - mi 

l / , when mi = m 3 & m-? = m 4 

V mi + m2 


m .4 


m 3 + m 4 
2 m 2 

mi + m 2 
m 4 

m 3 + m .4 

2 m 2 

mi + m 2 


1 ± 


m 2 m 3 


IT 


m2in 3 


(172) 


(173) 


(174) 


mi m. 4 , 

or 0 , when mi = m 3 & m 2 = m 4 , 


(175) 


mi 771-4, 

or 0 , when mi = m 3 & m 2 = m .4 , 


(176) 


and 


|p 4 (0i4 = tt/ 2)| = ±|pi| 


m .4 


m 2 m 3 


- 1 


m 3 T m .4 V mim 4 
= 0 , when mi = m 3 & m 2 = m .4 . (177) 

Note that solutions which give negative |p| are unphysical and must be discarded. 

(e) When mi = m 3 and m 2 = m. 4 , we found the physical solutions were 


|P3(6 i 13 = tt/ 2)| = |pi| 


I m 2 — m-i 
mi T m -2 ’ 


(178) 
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and 


|p 4 ( 0 i 4 = ?r/2) | = 0. 


(179) 


This says that the target can never be scattered into 9 i 4 = n/2. The projectile can only 
be scattered into 6 13 = n/2 if m 2 > rn \ . This makes sense. It is saying that the projectile 
can only be scattered into #13 = 7r/2 if it is lighter than the target. 

(f) When mi = 7773 and m 2 = 7774 , we obtained 


and 


|P3(013 = 0)| 
|P3(0i3 = T)| 
| P4 (6»14 = 0)1 


Pi 


or 


iPil or 

2777-2 


Pi 


7771 + 7772 


Pi I 
|Pl 


777-1 — 7772 
777-1 + 7772 
7772 ~ 7771 
7771 + 7772 


or 0, 


| P4 (6 , 14 = 7T ) | 


Pi 


27772 

7771 + 777-2 


or 0. 


(180) 

(181) 

(182) 


(183) 


We can never have the magnitude of a vector being a negative quantity, and so these 
solutions must be dismissed as unphysical. Also, we have already decided above that the 
cases |p 3 (#i 3 = 0)| = |pi| and |p 4 | = 0 correspond to the miss solution. Discarding the 
unphysical and miss solutions, we are left with 


and 


|P3(013 = 0)| 
|P3(013 = 7r) I 


Pll 

777-1 — 7772 

for 777-1 > 7772 , 

(184) 

777-1 + 777-2 ’ 
777-2 “ 777-1 

Pll 

for 777-2 > 7771 , 

(185) 

777-1 + 777-2 ’ 


|p 4 (014 = 0)1 


Pi 


27772 

7771 + 777-2 


(186) 


where we have indicated in parentheses the allowed mass values for physically acceptable 
solutions. Writing |p* | = 777 j|vj| gives (with rn 1 = 7773 and 777-2 = tt 7 4 ) 


|V 3 (013 

= 0)1 = 

l v ll 

7771 — 7772 

for 7771 > 7772 

7771 + 7772 ’ 

7772 — 7771 

|v 3 (6 , 13 

= 7T)| = 

l v ll 

for 7772 > 7771 

7771 + 7772 ’ 

|V 4 (01 4 

= 0)| = 

V 1 

27771 


7771 + 7772 



(187) 

(188) 

(189) 
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For mi » m2, we have 


(190) 

(191) 


|v 3 (013 = 0)1 « |vi| 

|v 4 (014 = 0)1 « 2|V! 

That is, the projectile continues on at the same speed and the target moves off at twice 
the speed of the projectile, both moving in the same direction as the original projectile. 

For m i = m2, we have 

|v 3 (0i 3 = 0 or 7r)| « 0, (192) 

|v 4 (0i4 = 0)| « |vi| . (193) 

That is, the projectile stops and the target moves off at the speed of the projectile, in the 
same direction as the projectile. 

For mi « m 2, we have 


|v 3 (0i 3 = vr)| « |vi| , (194) 

|v 4 (0i 4 = O)| « 0. (195) 

That is, the projectile bounces off from the target with the same speed but in the opposite 
direction, and the massive target remains stationary. 


3.2.3 Qualitative understanding of the |p 3 j| solutions 


We now develop a qualitative understanding of the two solutions, namely for a fixed angle 0i 3 j, 
there are two solutions for |p 3 ^ | . Consider the case of a billiard ball collision where a projectile 
billiard ball scatters off a target ball. In the above example, this is the case where mi = m3 and 
m2 = m 4 (elastic billiard ball collision), which was part (b) of the example. For that case, the 
momentum of particle 3 reduced to 


Write this as 


|P3Z| = |PlZ| 


mi cos 0i3 1 ± 


- sin 2 013/ 


mi + m2 


IPial 

IPszl 


= Pii 


= Pii 


mi ( cos 013/ + y — I - sin- 0i 3 / 
mi + m2 


mi cos 013/ - J - sin 2 0 i3 / 


mi + m2 


(196) 


(197) 

(198) 
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The solution |p 3 J allows for a negative value. This solution is obviously unphysical and may be 
discarded. For |p 3; | to be positive, we must have 


COS 0 i3 * 

> 

,H-sin 2 0 13 * 
y mf 

5 

cos 2 013 * 

> 

m l ■ 2 a 

— 2 - sin 0 i 3 * = 
mf 

2 

mf 

+ COS 2 013* - 1 
mf 

=» 1 

> 

2 > 
mf 


=>■ m\ 

> 

m 2 . 



Only if this is true, do we get both solutions |p^| and |p^|. Otherwise, we only get the one 
solution |p^|. So given that mi > m 2 , the question remains as to why we get two solutions |p 3 J 
at the same angle $ 13 ; for a fixed incident energy. Note that the two solutions are distinguished 
by their outgoing energy, or speed, or momentum magnitude. This is illustrated in figure 5, 
where a more massive projectile scatters off a less massive target. It can clearly be seen that 
the same scattering angle can occur for different impact parameters. If the projectile hits the 
target at small impact parameter then it will impart a relatively large amount of its momentum 
to the target and emerge with a relatively small speed at angle 9. 

On the other hand, if the projectile hits the target at a large impact parameter then it will 
impart a relatively small amount of its momentum to the target and emerge with a relatively 
higher speed at the same angle 9. Thus, the reason as to why there are two different momenta 
associated with the same scattering angle is because there can be two different impact parameters 
associated with that angle. Obviously, if we do an experiment and actually roll one billiard ball 
into another then there will be only one emerging momentum depending on the collision impact 
parameter. However, when we calculate a cross section, we must sum over all impact parameters. 
Thus, a cross section calculation will be the result of many billiard ball collisions at a variety 
of integrated impact parameters. The fact that the same lab angle gives rise to two different 
lab energies is also related to transformation from the cm frame to the lab frame. Later, we 
shall find that two different angles in the cm frame can give rise to the same angle in the lab 
frame. From figure 5, it is clear that the two different impact parameters will correspond to two 
different angles in the cm frame. 

3.3 Relations between initial and final momenta in lab and cm frames 

In this section, we develop some further useful results, involving relations between the initial 
and final momenta in the lab and cm frames as well as relations between angles. 

3.3.1 Lab frame 

Consider the elastic reaction 


1 + 2 — > 3 + 4 , 
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Figure 5: Elastic billiard ball collision where the mass of the projectile is greater than the mass 
of the target, b is the impact parameter. The balls are imagined to have the same density and 
so the more massive ball is drawn larger. (If the balls are of different density and have the same 
size, the dynamics shown in this figure are not altered. 


where the target particle 2 is at rest in the lab and the projectile particle 1 is fired into the 
stationary target. We will derive a formula relating the momenta of particles 3 and 4 with 
the relative angle $34 between particles 3 and 4. The formula will also contain the masses. 
This is much easier to work out by squaring the conservation of momentum equation, rather 
than breaking it up into components. After deriving this formula, we will show that it predicts 
#34 = 90° when all the particle masses are equal. 

We proceed as follows. Conservation of momentum is 

P 1 + P 2 = P3 + P4- (199) 


Squaring gives the non-relativistic Mandelstam variable, 

s = (P 1 + P 2) 2 = (P3 + P4) 2 

= I P3 1 2 T | P4 1 2 T 2 1 p 3 1 1 p 4 1 COS 0 34 
= |pu | 2 , in the lab frame where P2 1 = 0 . 

Conservation of kinetic energy can be written (with T 21 = 0), 

m 1 2 m 1 2 |_ |2 

|Pl[| _ | P3 1 I P4 1 

2mi 2 777-3 27774 

so that the Mandelstam variable becomes 

( I P3 1 2 . |P4| 2 \ 1 ,2 , I |2 , 01 II I a 

mi\ 1 = P3 + P4 + 2 P3 P4 COS 034 , 

\ 7774 J 


( 200 ) 


(201) 


(202) 
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which is the desired formula. When m\ = m 2 = m 3 = 777 . 4 , this becomes 

0 = cos 9 U =► 6> 3 4 = 90° . (203) 

3.3.2 Center of momentum frame 

Again, consider the elastic reaction 

1 + 2 — > 3 + 4 

but now we will analyze the reaction in the cm frame. Assume that we know the value of the 
momentum of particle 1 in the cm frame to be pi c . We can obtain some useful results. 

(a) We will work out a formula for |p3 C | as a function of |pi c | and show that it does not depend 
on any angle. 

(b) We will do the same for |p 4 C |. 

(c) We will show that the scattering angles are related by 0i 3c = 024c- 
We proceed as follows. 

(a) The non-relativistic Mandelstam variables are 



s = (P 1 + P 2) 2 = (P3 + P4) 2 , 

(204) 


t = (pi - P 3) 2 = (P4 - P 2) 2 , 

(205) 


u = (pi - P 4) 2 = (P3 - P 2) 2 , 

(206) 

and conservation of kinetic energy is 



2 2 2 2 

Pi _j_ P 2 _ P 3 Pi 

777-1 ^2 m 3 777-4 

(207) 

Now, specialize to the cm frame, where 



Pic + P2c = 0 = p 3c + p 4c , 

(208) 


|pic| = |P2c| = | Pic | and |p 3c | = |p 4c | = |p/ c | . 

(209) 

Eliminate p 2c and 

P 4 C | from the conservation of energy equation, giving 



Plc(— + — ) = P3c(— + — ). 

m 1 777-2 m 3 777-4 

(210) 


2 ,777! + 777-2 2 777 3 + m 4 

Plc( ) = P3c( ) > 

mi 777-2 7773 7774 

(211) 


„2 „2 m 3 m4(mi + m 2 ) 

P3c ^ >lc 777i?772(m3 + m4) 

(212) 
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Define 


2 _ m3m4(mi + 777. 2 ) 
77717772(7773 + 7774) 


which gives 


|P3c| = ^|Plc| , 

or 

|p/c| = +|p*c| , 


(213) 


(214) 


(215) 


which is the desired result, giving the relation between the momentum of particle 3 in 
terms of the incident momentum. No angle appears in these formulas. 

(b) In the cm frame, |pi c | = |P2c| = |p»c| and |p 3 C | = |p 4 C | = |pj c | so that 

|P4c| = M|P2c| 

= MlPicI ■ (216) 

(c) The non-relativistic Mandelstam t variable is 

t = (pi -p 3 ) 2 = (P 4 -P 2) 2 
= Pi + p| - 2|pi||p 3 | cos 6 >i 3 

= P 2 + P 4 ~ 2|p 2 ||p4| cos 024 • (217) 

In the cm frame, this becomes 

t = Pi + P/ - 2|pi||p/| cos #13 
= Pi + P/ - 2|pi||p/| COS024 , 

=+ COS 013c = COS 0 ' 24 c , 

=+ #13c = #24 c • 


4 Relativistic two body kinematics 

If the projectile or center of momentum frame moves fast enough, then two different angles in the 
center of momentum or projectile frames get boosted to only a single angle in the lab frame. In 
the lab, these two boosts will be distinguished by different lab energies or momenta. This is very 
important for 3 - dimensional space radiation transport codes which require cross sections in the 
lab (spacecraft) frame, which must be added incoherently from the two center of momentum or 
projectile frame angles. For a 2 - body final state, the formula which gives this double - valued 
lab momentum in terms of the lab angle, is given in equation (272). Other results, such as 
relating Mandelstam variables to energy and momentum variables, are also derived. 
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4.1 Introduction 

In this section, we wish to consider reactions of the type 

1 + 2— >3 + 4 + 5 + -- - , 

where the numbers denote the particles and the initial state always consists of just two bodies, 
whereas the final state can contain any number of particles. A 2 - body final state reaction is 

1 + 2 — >■ 3 + 4 , (218) 

whereas a 3 - body final state reaction is 

l + 2-^3 + 4 + 5. (219) 

The reaction kinematics are the same as shown in the previous figures 3 and 4. We are interested 
in calculating angular and spectral distributions for any of the final state particles. Thus, our 
primary interest is in the energy or angle of the jth final state particle, 

Ej or Oj , 

which give us the cross sections of primary interest, 

da da 

or . 

dEj dOj 

If Ej and Oj are dependent on each other, then and ^ are also dependent on each other 
and the doubly differential cross section, dE a e . has no meaning. In this case, the relationship 
between Ej and Oj is used to convert between dE and ^ . 

4.1.1 3 - body final state 

When there are 3 or more particles in the final state, then Ej and Oj are independent of each 
other. One can confirm this by working out the kinematics. In that case, one can form the 
differential cross sections jff, fjf- and dE f . 

4.1.2 2 - body final state 

This situation is more complicated. In the cm frame, Ej and Oj are independent, but still one 
cannot form fjjf - , because in the cm frame, Ej is fixed by the initial energy and the masses; see 

eqs. 235 and 236. However, Oj is undetermined and so one can form fjf- in the cm frame. In 
the lab frame, Ej and Oj are functions of each other and so jff and fjff are equivalent to each 
other. The proof of these statements is the main topic of this section. 
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4.1.3 Elastic and inelastic scattering 

Elastic scattering is defined as that in which kinetic energy and mass are conserved [Griffiths 
92], For the reaction with a 2 - body final state, as in equation (218), this will mean that 


Ti+T 2 = T3 + T4, (220) 

m\ + m 2 = m 3 + 777,4 . ( 221 ) 

If particles 3 and 4 are different from particles 1 and 2, then it will be very unlikely to find 
particles with just the right masses so that 7773 + 7774 = 7774 + m 2 . Thus, elastic scattering will 
almost always mean that ([19] ( p. 31), and [20] (p. 403) 


7771 = rn 3 , 7772 = ™4, |Pic| = |p/ c | = |p| (222) 

Thus, for all intents and purposes, we have the following statement. 

Elastic scattering occurs when the particle identites do not change. 

Examples are 


N + N -4 N + N , 


and 


n + N tt -{- N . 

An example of an inelastic reaction is 

N + N-+ N + A . 

This can be described with a similar formalism, simply by allowing particle masses to change 


[7] (p. 100). 

4.2 Mandelstam variables 

For the reaction 

1 + 2 -> 3 

the Mandelstam variables are defined as 

s = (P1+P2) 
t = (p\ ~ P3) 
U = (pi- pa) 


{P 3 +Pa ) 2 , 

(223) 

{P 2 -Pa ) 2 , 

(224) 

to 

1 

CO 

to 

(225) 
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Note that these are Lorentz invariant. An important result for the Mandelstam variables is that 
the sum of the variables is the sum of the masses squared, 

s + t + it = m 2 + ml + m 3 + m| . (226) 

4.2.1 Use of the s variable 

In the cm frame, y/s is the total energy, and 

\fs = E\ c + E2c • (227) 

The magnitude of the cm momenta can be written as 

| Pic | = IPicI = |P2c| , (228) 

|p/c| = |P3c| = |P4c| • (229) 


Defining 


A jk = A (s, = (s — m 2 — m\) 2 — ^m 2 m\ , 
one can then prove the following important results. 


Pic 



(230) 


(231) 


I A 34 


E\ c = 


E 2c — 


E3C — 


E Ac — 


C \ — 

V 4s 

1 

S + 

2 

m\ — 

ml 


\/4s 


s + 

2 

m 2 - 

2 


a/4 s 


s + 

mj - 

2 

m% 


a/4 s 


s + 

2 

m\ - 


a/4 s 


(232) 

(233) 

(234) 

(235) 

(236) 


Note that the energies and momenta of the final state particles, 3 and 4, in the cm frame are 
fixed. That is, they depend only on the masses and incident energy y/s. They do not depend on 
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any angles. This means that one cannot form the spectral distribution in the cm frame, (j 
is the label for particles 3 or f.) Conversely, therefore, the scattering angles 9 13 and 024 cannot 
depend on the momenta or energies of the final state particles. 

4.2.2 Use of the t variable 

Conservation of 4-momentum is 

P1+P2 = P3+P4, (237) 

and so the t variable can be defined and expanded as 

t = {pi -p 3 ) 2 = {P2 -Pi) 2 

= m\ + m\ - 2 EiE 3 + 2|pi||p 3 | cos 0 i3 

= m 2 + m\ — 2U2-E4 + 2|p2 1 IP4I cos 024 • (238) 

Therefore, a t distribution, (§r is equivalent to an angular distribution, %. Because the angles 
013 and 024 appear directly in t, it is useful in obtaining the relationship between Ej and 0 ? as 
mentioned above. This is discussed below. 

4.2.3 Equal mass particles 

Now, consider the case when all the particle masses are equal, i.e. mi 
One can show that | p* c | = |p/ c |- Label 

k = | Pic | — | p/c| j 

and let the cm scattering angle between particles 1 and 3 be 

013c = 9 ■ 

The following results can be shown ([2] (p. 102). 

s = 4 {k 2 + m 2 ) , 

t = —2k 2 {l — cos 0) = — 44’ 2 sin 2 (0/2) , 

u = — 2k 2 {l + cos 0) = —4k 2 cos 2 (0/2) . 

4.3 Relations between Ej and 9j 
4.3.1 Center of momentum frame 

In the cm frame, we have 

pi + p 2 = 0 , and p 3 + p 4 = 0 , (244) 


= m 2 = m3 = mi = m. 

(239) 

(240) 

(241) 

(242) 

(243) 
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showing that the relative angle $34 between the final state particles 3 and 4 is $34 = 180°. Thus, 
from figure 4, we see that 


Also, 


013 c — 024 c = 0 • 


(245) 



|Pi| = IP2I = |Pi| and |p 3 | = |p 4 | = |p/| . 

(246) 

We can 

use the t variable to show that 



s + ml - m\ 

E 3 c- r— 

V4s 

(247) 

This can be proved as follows. Use of equation (238) gives 



t = (pi -P3 ) 2 = {P 2 -P4) 2 

= ml + — 2 E 1 E 3 + 2 \pi\\p f\ cos 9 

= — 2 E 2 E 4 + 2|pj||p/| cos9 , 

(248) 

and the 

cos 9 term cancels out on both sides, so that 



ml + rr?J — 2 E 1 E 3 = rn'l + m\ — 2 E 2 E 4 . 

(249) 

Using the results 



Ei + E -2 = \/4 s , 

(250) 

and 

s + — m| = mf + — m\ — m\ + 2 E 2 (E± + E 2 ) , 

(251) 


which gives equation (247). Thus, we have used the t variable to evaluate the relation between 
Ej and 9j in the cm frame. We have found that the relation is trivial, in that E 3 does not 
depend on $13. We use this same technique below to find the (non-trivial) relation between Ej 
and 9j in the lab frame. 

4.3.2 Lab frame 

Conservation of 4-momentum gives 


Pi +P2 

= P3 + Pi , 

(252) 

t = {Pi -P3) 2 

II 

to 

1 

to 

(253) 

t = ml + ml - 2 EiE 3 + 2|pi||p 3 | cos 0i 3 

= rn'l + ml - 2m 2 E4i , 

(254) 
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or 


t = m\ + 777.3 — 2_E ] £3 + 2y £ 2 — vn\\jE\ — m 2 cos $13 . (255) 

Eliminating E±i, using conservation of energy, £ 4 ; = E\ + m 2 — £3 in the lab frame gives 
equation(254) as [10] (p. 102) 

E 31 (Eh + m 2 ) - |pu||P 3 i| cos 6 »i 3 ; = Eum 2 + + m\ + m\ - m\) 

= ^(s + mi — 777 . 4 ) 5 (lab frame), (256) 

with | p 1 / 1 = yj. Ef t — ml and | P 37 1 = yj E^ — m 3 . We have also used s = m 2 + m 2 + 2m 2 Eu. This 
is in agreement with references [10] (equation 4.65) and [7] (p. 127). Recall that when we did 
this analysis for the cm case, the cos 9 term canceled out on both sides. However, here for the lab 
case, no such cancelation occurs so that cos 9 remains, and this will give rise to the non-trivial 
relation between Ej and 9j. We now solve this important quadratic equation. Defining 


a = ^(tt7 2 + ml + ml - m\ + 2m 2 En) = ^(s + mj - 7774 ) , (257) 

b = Eh + 7772 , (258) 

c = | Pul cos 0i3i = \J £ 2 Z — rn\ cos 0 m , (259) 

enables equation (256) to be re-written as 

bE 3l - c|p 3 /| = a , (260) 

which is written either in terms of £ 3 ; as 

bE 3l - c\J Eh - = a , (261) 

or in terms of |p 3 j| as 

b \/ \P3i\ 2 + m 3 ~ c|P 3 /| = a . (262) 

Firstly, we solve equation (261). Squaring gives 

( 6 2 — c 2 )E 2 t — 2abE 3 i + (a 2 + c 2 m 3 ) = 0 , (263) 

which has the solution 

<76 ± c J a 2 — m 3 (b 2 — c 2 ) 

E 3l = , (264 ) 
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or 


E-si — 


ab± Ipul cos Oizi 
62 


\Ja 2 — m 3 ( b 2 

- p 2 u COS 2 013 1 


p li cos 2 e l 3 i) 


(265) 


with |pu| = Ey — m 2 . This equation agrees with Byckling [14] (equation 74) and Jackson [20] 
(equation 12.53). E 3 i has maximum and minimum values. Based on energy conservation alone, 
upper and lower bounds are given by 


4r er = m 3 , (266) 

and from energy conservation 

E 3 = Ei + E 2 — E 4 , 

E 3 i = Eh + m2 — En , 

which implies that 

E 3 f pC1 = Eu + m 2 - En = Eu + m 2 - mi . (269) 

Note that if m .2 = mi, then E ^ ax = Eu, which makes sense. The actual maximum and minimum 
values of energy are determined by conservation of both energy and momentum, and are found 
be setting 0 13 ; = 0 in equation 265. Equations (256 - 265) express the fundamental relation 
connecting E 3 i to di 3 1 . Given the fact that such a relation exists clearly shows that the angle is 
determined by the energy and vice versa. Thus, one cannot form ^ and ^ independently. ^ 
and ^ are equivalent to each other. Note, however, that we were unable to form ^ in the cm 
frame because the final momenta are constant and are not related to the angle. Secondly, we 
solve equation (262). Squaring gives 

(c 2 - 6 2 )|p3z| 2 + 2ac|p 3 *| + (a 2 - 6 2 mJ) = 0 , (270) 

which has the solution 

—ac ± bJa 2 + m?(c 2 — 6 2 ) 

|P3i| = , (271) 


(267) 

(268) 


or [21] (p. 52) 


P3Z 


-a|pn| cos 0i 3 ; ± b^j a 2 + m\ [ |piz| 2 cos 2 0 13 ; - (Eu + m 2 ) 2 ] 


IPizI 2 cos 2 0i 3 ; - (Eu + m 2 ) 2 


(272) 
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Physical solutions are found more readily using the |p 3 /| solutions in (272) rather than the E 31 
solutions in (265). The reason is as follows. The quantity 


«3 c = 


Pci 

P3c 


(273) 


is the speed of the cm frame with respect to the lab divided by the speed of particle 3 in the 
cm frame. As pointed out by Jackson [20] (p. 403), two roots are allowed in equation (265) for 
ct3 C > 1 but only one root for a^ c < 1 . Also, a physical constraint is that the solutions must 
obey E'3i > m3. However, for certain cases one obtains both roots no matter what the value 
of «3 c, even though in both cases E 31 > m 3 . Clearly, E 31 > m 3 is not enough to constrain the 
physical solutions [10]. The physical constraint on the magnitude of momentum is |p 3 ;| > 0. 
In the cases cited above, when one obtains two E31 roots for ct3 C < 1 , one finds that equation 
(272) gives one of the two |p 3 /| solutions as negative. (Of course one obtains the E31 solutions 

from |p 3 ;| by substituting into E31 = |P3Z | 2 + m|, using both the positive and negative values 
of |p3; | .) The negative solution is clearly unphysical and therefore one only has one physical 
solution when «3 C < 1 . In other words, the use of the |p3/| solutions in (272) rather than the 
E31 solutions in (265) enables us to immediately eliminate the unphysical solutions and we then 
see that two roots are allowed for ct3 C > 1 but only one root for ct3 C < 1. 

An example of an unphysical case discussed above occurs with the input parameters m\ = 
m3 = l,m2 = 7714 = 2, En = 2mi,9\3i = 7 t/ 7 , resulting in ot3 C = .625, = 1.86344, E^ = 

1.08538, |p^| = 1.57238, Ipj^] = —0.42196. As a final point, note that equation (270) can be 
written as 


A\p3i\ 2 - 2B\p 3l \ -<7 = 0, 

which has the solution 

B ± \JB 2 + AC _ B±D 
|PM| - A = 

where 


D = 

VB 2 + AC , 


E = 

E11 + m 2 , 

mj + 2 m 2 Eu) = ^(s + m| 

K = 

-{ml + ml + ml- 

A = 

E 2 - |pn| 2 cos 2 6*13/ 

= s + Pn| 2 sin 2 9 \ 3 i , 

B = 

K\pn\ cos 6*13/ , 


C = 

K 2 - m\E 2 . 



C can be re-written as 

4 C = (s + m 3 - m l) 2 - 4m|(£'n + m 2 ) 2 


(274) 


(275) 


(276) 

(277) 

(278) 

(279) 

(280) 

(281) 
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= 

(s + 777§ - ml) 2 - 

TYlo , 9 2\ 2 

^(5 + 7772-777.4) , 
7772 

(282) 

or 

4(7 7772 = 

2 / 1 2 2\ 

?77 2 (S + 777 3 — 777 4 ) 

2 2/ 1 2 2\2 

— m 3 (s + m 2 — m 1 ) . 

(283) 

Now, it turns out that 

03c > 1 

iff C < 0 , 

=>2 solutions, 

(284) 

«3c < 1 

iff C > 0 , 

=> 1 solution. 

(285) 


Here, the mathematical symbol iff means “if and only if’. « 3 C is given in equation (445) as 


«3c 


A 

/%c 


s + — m\ 

S — 771? + 777? 



Note that C(s) = 0 has the solutions 


•si 

S2 


7772(771.4 - 777.3) + m 3 (ml ~ ml ) 
7772 + 7773 

7?72(777.4 — 777|) — 777.3(777,4 — m^) 
7772 - 777.3 


(286) 


(287) 

(288) 


and a 3 c(s) = 1 has the same solutions, showing that C = 0 iff « 3 C = 1, which is consistent 
with (285). One can now vary s and study the behavior of C and a^ c . One can check (285) 
numerically by doing a parametric plot of C(s) and a 3 C (s), letting s vary. 


4.3.3 Proof of equation (285) 

First, note the following. 


A 

> 

0 

always, 

(289) 

B 2 

> 

0 

always, 

(290) 

I< 

> 

0 

always. 

(291) 


Equation (291) can be seen by noting that the minimum value of K = ^(s + 7773 — 777 . 4 ) occurs 

at the minimum value of s, which is at threshold where ^threshold = (m 3 + m^) 2 . (Note that 

s is Lorentz invariant so that this relation is true in any frame.) Thus, the minimum value 
of K is /^threshold = ^[( 777-3 + ?774) 2 + m 2 — 777 2 ] = 7773(7773 + 7774 ), which is always > 0. Given 
D = V B 2 + AC, it follows that (with \B\ = Vi} 2 ) 

D > \B\ iff C > 0, (292) 

D < \B\ iff C < 0. (293) 
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Secondly, note that we need to consider forward and backward scattering separately because of 
the following. 

Forward scattering: 

-|<0i3 i<\ => cos e 13i > 0 => B > 0 . 

Backward scattering: 

#13/ or 013/ > | => cos(9i 3Z <0 => B < 0. 

Now, write the solutions as 


\ -)- \ B + D 

Ips/I = A ? 

(294) 

B — D 

(295) 

Ip 3 /I = A > 


and note that 


D > 0 always. 

i) First, consider the case of forward scattering, i.e. B > 0. 

B > 0, A > 0 always, D > 0 always, => |p^| > 0, 

and, with B > 0, A > 0 always, equations (292) and (293) imply that 

|pjj| > 0 (allowed physically) if |B | > D => C < 0, 

|pgj < 0 (not allowed physically) if D > |B| => C > 0. 

Thus, 

if C < 0 => both roots |p^)| are allowed, 
if C > 0 => only one root |p^| is allowed. 

ii) Now, consider the case of backward scattering, i.e. B < 0. 

B < 0, A > 0 always, D > 0 always, => |p^| < 0 

and, with B < 0, A > 0 always, equations (292) and (293) imply that 

IPgjl > 0 (allowed physically) if D > |B| => C > 0, 

IP 3 /I < 0 ( n °t allowed physically) if D < |B| => C < 0. 
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Thus. 


If C < 0 => both roots |p^| are not allowed. 

If C > 0 => only one root |p(j)| is allowed. 

iii) Now, combine the results for both forward and backward scattering and summarize as follows. 

If C < 0 => both roots | p^| are either both allowed (forward scattering), 
or both not allowed (backward scattering) . 

If C > 0 => only one root |pj| is allowed, 

(for both forward and backward scattering) . 


4.4 Further results 

4.4.1 Momenta in the cm frame 

We wish to prove that |py c | = J ^ and |pj c | = J Note that using 


F 


any frame 


12 , 


(296) 


these equations imply that 


F c = \fs | Pic | = (Ei + E 2 ) |pic 


(297) 


We proceed with the proofs as follows. Eliminate p 4 in terms of the incident energy. Write 
the Mandelstam variable 


s = (pi +P 2) 2 = (P3 +Pi)' 2 

= (E\ + E2) 2 — (pi + P2) 2 = (E3 + E4) 2 — (P3 + P4) 2 . ( 298 ) 

In the cm frame, Pi + P2 = 0 = P3 + P4 and |pi| = | P 2 1 = | p* c | and | P3 1 = |p4 1 = |p/ c |- Thus, 

s = (Ei + E 2) 2 = (£3 + E 4) 2 = £3 + E 2 + 2 E 3 E 4 

= P/c + m \ + Pfc + ml + 2^j (p 2 c + ml)(p 2 fc + mf) , (299) 

and 

s ~~ m 3 — m \ ~ 2p/ c = 2 \/(P/c + m 3)(P/c + m 4) • ( 30 °) 

Squaring both sides gives 

(s - ml - ml) 2 - 4p 2 fc (s - m| - m\ ) + 4p} c 

= 4[p j c + P}c( m 3 + m l) + mlml] , (301) 
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and 


which gives [10] (p. 100) 


Similarly, 


4p 2 c s = (s — 777.3 — 777 2 ) 2 — 4m 2 777. 2 
= A(s, ml, ml) = A34 , 


P fc 



Pic 



(302) 


(303) 


(304) 


4.4.2 Energies in the cm frame 

Energies E\ c , E^ c and E 4 C may be expressed in terms of the Mandelstam invariant s, 



s = {pi+ P 2) 2 = {Ei + E 2 ) 2 - (pi + p 2 ) 2 . 

(305) 

In the cm 

frame (with pi = P 2 = p|), 



\fs = Ei + E -2 with E 2 = p 2 + 777 2 and Ef = p 2 + m 2 ■ 

(306) 

Thus, 

E 2 = E\ - mf + , 

(307) 

giving 

\fs — Ei = E-2 



= s]e\ - mj + 777-2 

(308) 


(a/s -Ei) 2 = Ef-ml + ml 

= s — \/4sEi + E 2 . 

(309) 

Thus, 

s — \ZAsE\ = — 7U 2 + ml , 

(310) 

giving 

S + 7T7- 2 — 77T, 2 
Elc - r— 

V4s 

(311) 
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The result for E 3c and E/± c are derived similarly [12] (p. 15). One can now calculate 


|pic| = \jE\ c -m\ , 

and the result agrees with the previously derived equation, 


I Pic = 


I Al2 
4s 


(312) 


(313) 


4.4.3 Bounds on Mandelstam variable t 

We wish to obtain formulas for to and t n . Specifically, we will prove that 
t 0 (U) = (. E lc - E 3c ) 2 - (|pic| T |P3 c|) 2 


rnf — m2 — m 3 + mf 


y/A s 


~ (|Plc| T |P3c 


1 


= —[(m? -ml -ml + mlf - (v^T v^) 2 ] • 


(314) 


The notation =p means that the equation for to has the — sign and the equation for t^ has the 
+ sign. We proceed with the proof as follows. The Mandelstam variable is 


(315) 


and we use 


t = (pj - p 3 ) 2 

= (E 1 - E 3 ) 2 - ( Pl - p 3 ) 2 


(P1-P3) 2 = |Pl| 2 + |P3| 2 - 2pi-p 3 

= |pi| 2 + IP 3 I 2 - 2|pi| |p 3 | cos 6» 

= |pi| 2 + I P3 1 2 - 2 |pi | |p 3 1 (1 - 2sin 2 °-) 

Q 

= |Pi| 2 + I P3 1 2 - 2 |pi | |p 3 1 + 4 |pi | |p 3 1 sin 2 - 
= (|pi| - I P3 1 ) 2 + 4| Pi 1 1 p 3 1 sin 2 ^ . (316) 

With the definitions, to = t(6 c = 0) and t n = t{6 c = 7r), we get 

to(U) = (Eic ~ E 3c ) 2 - (|pic| t |P3c|) 2 (317) 

as desired. For the proof of the second equation above, use the previously derived results, 


Ei c — 


s + m 2 — m 2 


\[A s 


and E 3c = 


s + nr 2 — m 2 


i/4s 


(318) 
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to give 


which results in 

t 


(E lc - E 3c ) 2 


m\ — m 2 — m 2 + m\ 
y/4 ~s 


{E lc - E 3c ) 2 - (pic - p 3c ) 2 

(Elc ~ E 3c ) 2 - ( | p r | - |p 3 |) 2 — 4|p 1 ||p 3 | sin 2 j . 


(319) 


(320) 


With the definitions to = t(6 c = 0) and t n = t(0 c = ir), we obtain equation (314) as desired. 


4.4.4 Momentum and energy thresholds 

We now obtain lab energy and momentum thresholds. Consider the 2 - body reaction 

1 + 2— >3 + 4. (321) 

We wish to prove the following results. 

(a) If the lab momentum of particle 1 is pu, we will derive expressions for s = (pi + P 2) 2 = 
(P 3 + P 4,) 2 in terms of E lh T u and pu. 

(b) We will obtain the threshold value of s needed to produce the final state particles. 

(c) We will obtain the threshold value of Tu needed to produce the final state particles. 

We now proceed with the proofs as follows. 

(a) Start with 

S = (pi+P2) 2 v 

= (Ei + E 2) 2 — (pi + P 2) 2 • (322) 

In the lab frame, P 2 = 0 and therefore E 21 = m 2 . Thus, in the lab frame 
s = ( Eu + ^ 2) 2 - P? 

= E11 + Trr^En + m| — pj 

9 9 0 n / O O 

= rn 1 + m 2 + 2m2Eu = m x + m 2 + 2 m 2 \! p\i + mf 
= (mi + m 2) 2 + 2m 2 T li . (323) 

(b) s is invariant. In the cm frame, s = [E\ + E 2) 2 = (E 3 + E 4) 2 and so the threshold value is 

'^threshold = (n ^ 3 + m. 4 ) (324) 

when the final state particles are all produced at rest. 
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(c) Substituting into the equation above, we get 


Til threshold — 


(m 3 + ^ 4) 2 - (mi + 777 2 ) 2 
2777,2 


(325) 


4.4.5 Relative scattering angle in the cm frame 


Now consider the relative angle between the two final state particles. Consider the reaction 
1 + 2 — > 3 + 4. We can obtain some more useful results. 

(a) We show that in the center of momentum frame, the particles 3 and 4 scatter at a relative 
angle of 6^ c = 180°. We show this is true in both the non-relativistic and relativistic 
cases. 


(a) Now, consider the lab frame. We will show that 

Pi = p! + P 4 + 2|P3| |Pd| COS 0347 . 


(326) 


where $ 34 / is the relative scattering angle between particles 3 and 4 in lab frame, and p 
refers to the particle 3-momenta. We show this equation is true in both the non-relativistic 
and relativistic cases. 

(c) Now, consider part b) for the case of elastic scattering, where all particle masses are equal. 
In the non-relativistic case, we will show that the lab angle is $347 = 90°. (Note: In the 
non-relativistic case, elastic scattering is defined as conservation of kinetic energy.) 

(d) Again, consider part b) for the case of elastic scattering, where all particle masses are 
equal. This time consider the relativistic case. Now, the scattering angle between the final 
particles is not 90°, but is instead given by 


cos 6*347 = 


r 3 r 4 

I P3 1 1 P4 1 


T 3 (Ti-T 3 ) 


\Jt% + 2777T 3V / (T 1 - T 3) 2 + 2m(T 1 - T 3 ) 


(327) 


where T 3 and X 4 are the kinetic energies of particles 3 and 4 and |p 3 1 and |p 4 j are the 
magnitudes of their 3-momenta. We prove this result for the relativistic case. Also, we 
show how this reduces to the non-relativistic result $347 = 90° for small speeds. (Note: 
In the relativistic case, elastic scattering is defined as conservation of kinetic energy and 
mass. Small speed is equivalent to mass being much bigger than kinetic energy.) 

We now prove these results as follows. 

(a) The cm frame, in both the non-relativistic and relativistic cases, is defined by 


P3 + P4 = 0 , 


(328) 


giving 


P3 = -P4 


(329) 
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implying that the relative angle in the center of momentum frame is 


9 Uc = 180° . (330) 

Alternative method: 


P3 = -P4 • 


Define 


Thus, 


_ 2 '2 2 

P3 = P4 = P / • 


Therefore, 


(P3 + P4 ) 2 = 0 


p! + pi + 2|p 3 ||p 4 | COS034C 

2p } + 2p j COS 6*34 c . 


2 P / 

COS $34c 
$34c 


-2py COS 034 C 
-1 

180° , in the cm frame . 


(331) 


(332) 


(333) 


(334) 

(335) 

(336) 


This is true in both the non-relativistic and relativistic cases. No assumptions have been 
made concerning mass or whether or not the scattering is elastic. 

(b) The lab frame in both the non-relativistic and relativistic cases is defined by 


P 2 = 0 . (337) 

Conservation of 3-momentum becomes 

Pi + P 2 = P 3 + P 4 (338) 

Pi = P 3 + P 4 , lab frame . (339) 

=> Pi = P 3 T p 2 T 2 1 P 3 1 1 P 4 1 cos $ 34 / . (340) 


(c) In the non-relativistic case, elastic scattering is defined by conservation of kinetic 

T\ = T 3 + T 4 , lab frame. 

P? = P§ + P| 

2m 2 m 2m 

= > Pi = P3 + P4- 

Substituting into equation (340) gives cos 6 * 34 / = 0 or 0 34 ; = 90°. 


energy, 

(341) 

(342) 

(343) 
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(d) In the relativistic case, total energy is always conserved, 


Ei + E 2 — E 3 + E 4 , (344) 

Ti + mi + T 2 + m 2 = T 3 + m. 3 + T 4 + m 4 . (345) 

In relativistic elastic collisions, kinetic energy, rest energy and mass are all conserved. 
Thus, for relativistic elastic collisions 



Ti + T 2 = T 3 + T 4 , 
777.1 + m 2 = m 3 + 7774 . 

(346) 

(347) 

The second equation is of course true by our assumption of equal 
lab frame, 

mass particles. In the 


Ti = T - 3 + T 4 lab frame 

(348) 

which is the same as the non-relativistic case. However, now the 
energy are different, namely 

expressions for kinetic 

\J Pi + m — m = 

\J p 3 T 777. — 777. T \J p 4 T 777 — 777. , 

lab frame , (349) 

\Jp\ + m + 777 = 

= Ei + m = 

\Jv > 3 + m + \Jp\ + 777. 
e 3 + e 4 , 

(350) 


where we could have obtained this equation right away by setting E 2 = m in equation 
(344). This all illustrates how the conservation of energy and kinetic energy work out in 
the relativistic case, and we have written this out for comparison to the non-relativstic 
results. However, to actually get a result in terms of the angle # 34 /, it is easier to work 
with 4-vectors. Thus, 


0 P 1 +P 2 ) 2 = (P3 + P4) 2 , (351) 

rri'i + ini + 2(EiE 2 - p 4 • p 2 ) = rn 3 + m\ + 2 (E 3 E 4 - p 3 • p 4 ) . (352) 

Let mi = rn 2 = niz = m 4 = m. This implies elastic scattering, because mass is conserved 
and as shown above, it implies conservation of kinetic energy. This gives 


E\E 2 - pi • p 2 = E 3 E a - p 3 • p 4 . 


The lab frame is defined with p 2 = 0, giving 


mEi 
cos 0 34 i 


E 3 E 4 - |p 3 ||p 4 | cos 0 34 i 
E 3 E 4 — mEi 
| P3 1 1 P4 1 


(353) 


(354) 

(355) 


Conservation of energy gives Ei+E 2 = E 3 +E 4 which, in the lab frame, is Ei+m = E 3 +E 4 , 
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or Ei = E 3 + E 4 — m. Substitution gives 


cos 6 * 34 1 = 


T 3 T 4 


(356) 


|P3||P4| 

This is not the final answer, however, because Ti = T 3 + T 4 giving T 4 = T\ — T 3 . Using 
p 2 = E 2 — m 2 = (T + m ) 2 — m 2 = T 2 + 2 mT , (357) 

gives 


p| = T 3 2 + 2 mT 3 , (358) 


and 


Thus, 


P 4 — Ti + 2mT/i — {T\ — T 3) 2 + 2m{T\ — T 3 ) . 


COS 6>34Z = 


T 3 T 4 


r 3 (Ti-T 3 ) 


(359) 


(360) 


IPs 1 1 P4 1 yT| + 2mT 3 v /(T 1 - T 3 ) 2 + 2m(Ti - T 3 ) 

In the non- relativistic limit, energy is dominated by rest mass energy. In other words, 

m » T . (361) 

Thus, 


cos 9 3 4 l 


giving 


n m - t 3 ) 

v /2mT 3 v / (2m(Ti - T 3 ) 
y /T 3 (Ti - T 3 ) 

2 m 

0 , for m » T , 


0341 ~ 90 c 


(362) 


(363) 


in the non-relativistic limit. 


4.4.6 Equivalence of |pj c | and |p/ c | 

Consider the reaction 

1 + 2-13 + 4. 
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Consider the case when all the particle masses are equal, i.e. m\ = m 2 = m3 = = m. One 

can show that | p* c | = |p/ c |- (See the previous results). Label 

k = |pi c | = |p/| , (364) 

and let the cm scattering angle between particles 1 and 3 be 

0i3c = 0 • (365) 

We will show the following results. 

(a) s = 4 (k 2 + m 2 ) , (366) 

(b) t = — 2A: 2 (1 — cos 0) = — 4& 2 sin 2 (0/2) , (367) 

(c) u = — 2/c 2 (l + cos 0) = — 4/c 2 cos 2 (0/2) . (368) 


To show these results, we proceed as follows. 

(a) Begin with 

s = ('P 1 +P 2) 2 

= P\ + P 2 + 2 Pi -P 2 

= ml + + 2(EiE 2 — P 1 .P 2 ) • (369) 

In the cm frame, Pi + P2 = 0 giving P1.P2 = — | P 1 1 2 = — p 2 c = —k 2 . Thus 

s = m 2 + m 2 + 2[y ( k 2 + m 2 )(k 2 + m 2 ) + k 2 ] 

= 2 m 2 + 2 (k 2 + m 2 + k 2 ) with m± = m 2 = m 
= 4 (k 2 + m 2 ). (370) 

(b) The Mandelstam variable is 

t = ( Pl - P3 ) 2 

= pl+pi-2pi-P3 
= 7nj + ml-2(E 1 E 2 -pi.p 3 ) 

= m\ + (pf + m 2 )(p 2 + mf) + 2|pi||p 3 | cos0i 3 . (371) 

In the cm frame, when m± = m> 2 = m 3 = 1/14 = m, we have |pi c | = |p 3c | = k, which gives 

t = 2m 2 — 2(k 2 + m 2 ) + 2k 2 cos 0 

= — 2fc 2 (l — cos0) = — 4L 2 sin 2 (0/2) . (372) 
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(c) Changing 3 — > 4 in the above gives 

u = (pi- Pa ) 2 

= — 2k 2 (1 — cos 0i4) 

= —2k 2 [l — cos(180 — 0 13 )] 

= — 2k 2 (1 + cos 0i3) 

= — 2k 2 (l T cos 9) = — 4/c 2 cos 2 (#/2) where 9 = 613 . (373) 

The relation 9 14 = 180 — $13 can be seen in figure 6. 
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Figure 6: The reaction, 1 + 2 — > 3 + 4, viewed in the cm frame. 


4.4.7 Relationship of s and pu 

For the reaction, 


1 + 2 — > anything , 


we now show that 


s = m\ + m 2 + 2rn2E\i 

= m\ + m 2 + 2m 2 \Jph + m 2 (374) 

where pn is the momentum of particle 1 in the lab frame. This is proved as follows. The 
Mandelstam variables is 


s = (pi +P2) 2 = pj + pl + %Pi-P2 
= to 2 + ml + 2(EiE 2 - P 1 .P 2 ) ■ 


(375) 
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The lab (or target) frame is defined as the frame in which the target is at rest, namely p 2 1 = 0. 
Thus, equation (374) is obtained. 

4.4.8 Threshold of s 

For the reaction, 


1 + 2— >3 + 4 + 5, 


we now derive an expression for the threshold value of s in terms of to, 3 , TO 4 , rri^ only. We 
proceed as follows. Our formula for the threshold kinetic energy was given in equation (325) 
and generalizes to 3 final state particles as 

_ (m . 3 + TO4 + to 5 ) 2 - (toi + m2) 2 

1 U = ^ 

2 m 2 

In the previous result, we showed that 



s = mf + m 2 + 2m 2 E n = m\ + m 2 + 2 m 2 {Tu + toi) 

= (toi + to 2 ) 2 + 2m2+u = (mi + m 2 ) 2 + (to 3 + to 4 + to, 5 ) 2 - (mi + to 2 ) 2 
= (to 3 + to 4 + to, 5 ) 2 . (377) 


An alternative derivation is 


s = ('P1+P2) 2 

= (p 3 + pi + P 5) 2 , (due to 4 — momentum conservation) 

= (to 3 + ?n 4 + to , 5 ) 2 , (378) 

because all final state particles are produced at rest in the cm frame at threshold. 


5 Lorentz transformations 

Cross sections are usually calculated in the center of momentum or projectile frames, but space 
radiation transport codes require cross section formulas in the lab (spacecraft) frame. One 
therefore must transform cross sections, angle, energies and momenta into the lab frame, by 
using Lorentz transformations. The general formulas for doing this are derived in this section. 

Suppose a frame S' is moving at velocity v relative to a stationary frame S, as shown in 
figure 7. We will always choose the longitudinal z direction to be in the direction of motion. 
The transverse directions x and y will always be perpendicular to the velocity. The Lorentz 
tranformations for transforming spacetime coordinates from one frame to another as shown in 
figure 7 are 


t' 


7 

- 7/3 


7/3 

7 


t 

z 



(379) 
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and 


x = x y = y , (380) 

where 

P = v. (381) 



Figure 7: 
Frame S moves relative to 


Reference frame S' moves at speed v relative to S. 

S' at speed —v, as shown in figure 8. The inverse transformation is 




(382) 


Thus, the inverse transformations are obtained by swapping primes and unprimes and changing 
the sign ofv. 


5.1 Lorentz transformations of energy and momentum 

Consider a frame S' moving at speed v relative to frame S. Now, imagine that a particle is 
moving at velocity u in frame S'. (See figure 2-4 from Tipler [22].) If the particle has energy 
and momentum E' and p' as measured in frame S', the Lorentz transformation gives the energy 
and momentum E and p measured in frame S. 

The coordinates t and x form a 4-vector (t, x) with an invariant squared length given by 
s 2 = t 2 — x 2 . Similarly, for the energy-momentum 4- vector (E, p), the invariant squared ’’length” 
is E 2 — p 2 = m 2 . The E and p obey Lorentz transformations identical to t and x: 


E' 


( 7 ~lP \ ( E \ 

V -lP 7 ) \Pz ) ’ 


(383) 


Px = Px , 


Py = Py ■ 
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(384) 




Figure 8: Frame S viewed from S'. 


It is often very useful to write this in terms of parallel momenta p\\ = p z and transverse momenta 
PTx = Px or pT y = Py We often just write px to denote either px x or px y ■ The direction parallel 
or transverse is defined in relation to the velocity of the frame S'. Thus, we write 


E' 

P\\ 


( 1 - 7/5 

V -7/5 7 



Pt = PT ■ 


(385) 


5.2 Transformation between cm or projectile frame and lab (target) frame 

Suppose we have quantities in the cm or projectile frames and we wish to transform to the 
lab frame. The cm frame moves at speed ( 3 c i relative to the lab frame. The corresponding 7 
factor is labelled as 7 c /. The projectile frame moves at speed ( 5 p i relative to the lab frame. The 
corresponding 7 factor is labelled as 7 p i. The Lorentz transformations are 


( E * ) = I 

V p \\* J 

and inverse transformations are 

(») - 
V P\\i J 

where 


1*1 

—7*1 ( 3 *i 


~l*if3*i W Ei 

l*i ) { P\\i 


, PT* = PTl j 


1*1 l*i(3*i \ ( E* 

l*i(3*i l*i j l P\\* 


, PTl = PT* , 


P\\=Pz = |p| cos 0 , 

Pt = |p| sin 6 . 


(386) 


(387) 


(388) 

(389) 


With the above notation, both cm and projectile frames are included. The notation is as follows. 
A quantity x* is the value of the quantity x evaluated in that particluar frame with 


* = c or * = p , 


(390) 


56 



and (3d is the speed of that frame with respect to the lab with 


(3*i = (3 c i or (3*i = (3 p i . 


(391) 


5.2.1 Evaluation of (3d and y d 

Consider the non-relativistic case. For 2 particles, the position of the center of mass frame, 
relative to an arbitrary origin, is defined via 


R = 


rri\Y\ + m 2 r 2 
m\ + m 2 


giving the velocity of the cm frame as 


(392) 


V = 


m ] v i + miv 2 pi + p 2 


mi + m 2 mi + m 2 

The velocity of the cm frame relative to the lab frame is obtained by setting p 2 = 0 to give 

Pi i 


(393) 


V d = 


mi + m 2 


(394) 


where pn is the momentum of particle 1 (the projectile) in the lab frame. The relativistic case 
follows similarly. For a single particle of mass m, p = ymv and E = 7 m giving 



= v 

(395) 


E 

1 = — , 
m 

(396) 


0 , £ 

II 

> 

(397) 

For a system of particles (Byckling [14] (p. 21) of total energy E = E\ + E 2 + • • 
momentum p = pi + p 2 + • • •, the invariant mass is 

• and total 

M 

III 

to 

1 

to 

II 

(398) 

which, for 2 particles, is [1] 



M =y/(E 1 + E 2 ) 2 - (pi + p 2 ) 2 = y/s . 

(399) 

The velocity of the cm frame is 



(3 = v 

P _ Pi + P2 H 

E E\ + E 2 + • • • 

(400) 

7 

E Ei + E 2 + • • • 
M ~ 

(401) 

1(3 

P Pi + P2 H 

M ~ J~s 

(402) 
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For a system of 2 particles, the velocity of the cm frame relative to the lab frame is obtained by 
setting p 2 = 0, to give [1] 


Pci = 

Pi/ 

(403) 

Eu + rn 2 

7 Cl = 

En + m 2 

(404) 

V~s ’ 

7 cl Pci 

Pi/ 

(405) 


with 

\fs = \J m\ + m| + 2 m 2 Eu , (406) 

where pn is the momentum of particle 1 (the projectile) in the lab frame. These results can 
also be written in terms of invariants. It is straightforward to show that 

yA^ = 2m 2 | Pl/ | , (407) 

and because A is invariant, it does not need a frame label, i.e. 


iPiil = 

Substituting into the above expressions gives 

Pd = 

Id = 


\Al2 
2 m 2 

([14] pps. 26 and 73) 
n/Ai 2 

S — mf + 1712 1 

s — mf + m| 

2m 2 y / s 


' yd Pci 



5.2.2 Evaluation of P p i and 

Clearly, 


Ppi 


'Ipi 


7 pi Ppi 


Pi/ 

Eu 

Eu 

m\ 
Pi i 
rri\ 


1 + 


Tu 

mi 


(408) 

(409) 

(410) 

(411) 


(412) 

(413) 

(414) 
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where T is the kinetic energy. Again, these results can also be written in terms of invariants. It 
is straightforward to show that 


Ppl 

1 pi 

7 pi ft pi 


\Al2 

2 2 ’ 
S — m\ — 777-2 

s — 


2mim2 



2mim2 


(415) 

(416) 

(417) 


5.3 Energy transformations 

The Lorentz transformation for energy from the lab (l) frame to the starred (*) frame is 


E * = i*i(E t - (3*ip\\i) 

= i*i(E t - fi*i\pi\ cos 6i) 

= l*i [ei - (3*i \] Ef - m 2 cos B^j . (418) 

The derivative dE*/dE\ is found by expressing cosOi in (254) in terms of Ei, and substituting for 
cosOi above to eliminate 9i and have E * solely as a function of £). The inverse transformation is 

Ei = i*i (e* + (3*i \J El -m? cos d^j . (419) 


5.4 Angle transformations 

The angle is obtained from 

Pz 


Thus, the angle of particle j is 

tan 6ji = ^ 
Pzjl 


PTj* 

l*lP*lEj* T i*iPz j* 

|Pj*l sin^j* 

l*i ((3*i Ej* + |pj*| cos Oj*) 


(420) 


(421) 


Defining a,j * as the speed of the (cm or projectile) frame relative to the lab divided by the speed 
of particle j in the (cm or projectile) frame 


a j* 


(3*1 

(3j* 


(422) 
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and using 


fij* 


Pi* 



we obtain 


tan Oji = 


sin Oj * 

5*1 (cos Oj* T (Xj*') 


(423) 


(424) 


See references [20] (p. 402), [14] (p. 42), [12] (p. 17), [19] (p. 26). The above equation is a 
complicated function of 6 because in general 


— oq ( Ej * ) — oy * (0 j*j , 


(425) 


To make this explicit, the equation is written 


tan Oji = 


sin Oj * 

7*/ [cos Oj * T otj*(0j*)\ 


(426) 


That is, in general, oij * is a function of Oj * , making tan Oji a complicated function of Oj*. How- 
ever, for the cm frame Ej c is not a function of 0j c , meaning that ctj c is not a function of 0j c 
[If] (pp. f2, 58). Also, for a 3 - body final state E JC is not a function of 0j c meaning that oy c is 
not a function of 0j c . For * = c, this is plotted in references [14] (p. 43), [20] (p. 403), [23]. For 
ctj * > 1, the function is double valued. That is, two different angles in the (cm or projectile) 
frame can give rise to the same angle in the lab frame for ay * > 1. However, the two particles 
can be distinguished by their energies, labelled for the cm frame as E± c [20] (p. 402). Later we 
will need to evaluate ■ Using sec 2 0 — tan 2 0 = 1 allows equation (426) to be re-written as 


cos Oji 


5*1 [cos 0j* 4~ dj * (0j * ) ] 



+ aj*(0j*)} 2 + sin 2 0j * 


(427) 


which is written in terms of 0j*, to be contrasted later with (433) written in terms of Oji. In 
evaluating , one needs to be very careful, because, as stated above, atj * is a function of 0j*. 

That is, in general, a.j * is a function of 0j* , making the evaluation of difficult. However, 

for the cm frame, Ej c is not a function of 0j c , meaning that ctj c is not a function of 0j c ([If] 
pps. f2 and 58). This is also true for a 3 - body final state. Thus, for the cm frame we have 
[14] (p. 59) 


d cos Oji _ 7*z(l + dj c cos 0j C ) 

d cos 0j c [ihiotjc + cos 0j c ) 2 + sin 2 0 jc \ 3 / 2 ’ 


(428) 


which is written in terms of 0j c , to be contrasted later with [14] (p. 59) written in terms of Oji . 
Using exactly the same technique that was used to derive equation (426), we may also derive 
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the following [14] (p. 42) 


tan Of = 

J* 


sin dji 
7*/ [cos Oji - af^Oji)] 


where 


% 


a t = 


Again, we have been careful to note that is a function of Oji, i 


zf, 

P*i 


i.e. 


a fi = 


(429) 


(430) 

(431) 


(432) 


This time the lab momentum and energy are double valued functions of the lab angle, as dis- 
cussed previously. This is not the case in the cm frame. There, the energy and momentum are 
independent of the angle. Now, equation (429) is a complicated function of 9ji because in general 
a.ji = aji(9ji). An alternative to equation (429) is to invert equation (426) directly which leads 
to [14] (p. 43, equation 8.32) 


cos 6j * 


— aj*(0j*)7*i tan 2 Oji ± \[T) 
1 + 7 h tan 2 Oji 


(433) 


which is written in terms of dji, in contrast with equation (427) written in terms of 0-,*. The 
term D is given by 


D = 1 + 7*; tan 2 9 j{ [ 1 - a|*(0j*)] , 


(434) 


and using [14] (p. 43) 


1 + 7 h tan 2 dji 


ill 


cos^ 


Vjl 


-(1 - 0h cos 2 Oji) 


(435) 


gives 


COS 


—Oij*{0j*)^h[l — cos 2 Oji } ± cos 2 OjiVD 

iU 1 - Ph cos2 °ji ) 


(436) 


This is still a complicated function of 0j* through ay* = ay* (61/*). However, it has a big advantage 
for the cm frame because, as discussed above, oy c is not a function of 0 ]C . The derivative is [14] 
(p. 59) 


d cos 0j c cos Oji ( a.j c ±\fD ) 2 

d cos Oji 7 2 ;(1 - /3 2 ; cos 2 Oji) 2 (±\/D) 
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2 


= (MY 1 

\|Pjc| / {±\fD cos 6 ji) 

_i_ I n |2 

= IPj;| , (437) 

r 7d\Pjc\{\Pji\ - EjiP d cos 6 ji) ’ 

which is written in terms of dji, in contrast with equation (428) written in terms of 0j c . Here, 
|Pj«| given by (272) and E jt = \J\Pji\ 2 + rrif where 

D = 1 + tan 2 Oji( 1 — Oj c ) . (438) 

5.4.1 Evaluation of a^ c and ol/^ c 

We have already obtained /3 c i in equation (409) above. Now, 

Pic = ^ • (439) 


Using the results from section 4.2.1, 


|P3c| = I P4c | = |P/c| = 



Ezc — 


s + m§ — m .4 


E4 c — 


s + m 2 — ru| 


which gives 


o _ 34 

P3c 1 2 2 > 

S + m.3 — m A 

a _ \/^34 

P4c — ,2 2 ) 

S + 777.4 — m 3 


which results in [14] (p. 73) 


G?4 c — 


S d - 777 2 m 4 / A 12 
S — 777 2 + 777.| V ^3 4 

S + ?77 2 - TTT . 2 M 12 

S — 777 2 + m 2 V ^3 4 


We see that a 3 C is the same as a 4 C except for the interchange 3 <-7 4. This gives exactly the 
same result as a in equation (12.51) of reference [20] and r in reference [19]. 
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5.5 Transformation of angular distributions 

Suppose we have in the cm or projectile frames and we wish to obtain ^yy in the lab frame. 
Clearly, 


da 

d cos Oji 


da d cos 6j * 
d cos Oj * d cos dji 


and using dfl = 27rd(cos 0) gives 

da d cos Oj * da 

dflji d cos Oji dQj* 


Using equation (428) gives 


da [7 *i( a jc + cos Ojc) 2 + sin 2 0j c ] 3 ^ 2 da 

dtlji 7 * 2 1 1 + a j c cos 0j C \ dQj c 


(447) 


(448) 


(449) 


which agrees with Joachain [19] (equation 2.104). This is written in terms of 9j c , in contrast to 
(450) written in terms of Oji- Equation (440) is one of our most important results. It tells us 
how to transform an angular distribution from the cm frame to the lab frame. The way to use it 
is as follows. Suppose you want the lab frame angular distribution -if- evaluated at a particular 
lab angle, say Oji = 45°. Then, evaluate the corresponding cm angle 0j c by substituting Oji = 45° 
into equation (433). This cm angle is simply substituted into the right hand side of equation 
(449), and the number obtained is the value of . Now of course, because of the double valued 

nature of (433), you might get two angles, say 0+ c and 0~ c . If this occurs, it simply means 

that there are two terms on the right hand side of (449) which must be added incoherently (i.e. 

ordinary simple addition with no interference terms). However, equation (449) is inconvenient 

in one sense. If we want the angular distribution -4ff- in the lab frame, then it would be very 

“jl 

nice to have it only as a function of the lab frame angle Oji. But equation (449) is written in 
terms of the cm angle 6j c . Using reference [14] (p. 59) gives 


da 

dTtji 


cosOji (« 3 C ± \/D) 2 da 

7 dC 1 - Pd cos e ji) 2 (±/D) dn jc 

( mV 1 da 

\ Ml ) (±V^COS 0 J7 ) dTljc 

M 2 da 

7dMI(M - EjiPdCosOji) dQ jc ‘ 


(450) 


This is written in terms of Oji, in contrast to equation (449) written in terms of 0 JC . However, 
the right hand side still contains which is a function of 0 JC . To eliminate this dependance, 

substitute equation (436) wherever cos 0j c appears in This will result in the right hand side 
of equation (450) being only a function of Oji. Of course, it will be double valued sometimes. 
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5.6 Double differential cross sections 


In order to transform an area element such as dxdy or a volume element such as dxdydz, we 
need to use the Jacobian defined as 


d(x,y) 
d(u, v) 


dx 

dx 

du 

dv 

dy 

dy 

du 

dv 


dx dy 
du dv 


dy dx 
du dv ’ 


(451) 


so that under a change of coordinates, an integral becomes 


J J f{x,y)dxdy = j J f[x(u, v),y(u, u)] V J dudv 


(452) 


or, in other words, 


dxdy = dudv . (453) 

d(u, v) 

The proof of this is usually given in advanced calculus books. If several transformations need 
to be carried out, we can use the chain rule properties of the Jacobian, namely 


d(x i ■■■x n ) _ d(x i • • • x n ) d(yi ■■■y n ) 
d{zi ■■■ z n ) d(yi ■ ■ ■ y n ) d(zi ■ ■ ■ z n ) 


(454) 


We are interested in relating the double differential cross section between the cm or projectile 
frames and the lab frame, as in 


d 2 a d 2 a d(Ej*,Qj*) 

dEjidQji dEj^dQj* d(Eji,Qji) 

which is the analog of equation (447). In the following analysis, 
are independent. This is true for 3 particles in the final state, but, as we saw previously, it is 
not true for only 2 particles in the final state. Re-write the Lorentz transformations (387) as 
(leaving off the index i) 


(455) 


it is essential that E[ and 9i 


or 


P* 

Ip* 


.E* 

II 

1 

P 

(456) 

P\\* 

= T*lip\\i - P*l E l) s 

(457) 

Pt* 

= PTl , 

(458) 


E * = 7 *;(£) — /3*z pz \cos9i) , 

(459) 

cos 9* = 7*K P« cos 9i - f3*iEi) , 

(460) 

sin 6b = |p;| sin 9i . 

(461) 
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Now our two independent variables are E and 6 (or Q). For the case of 3 or more particles in 
the final state, E is no longer a function of 6. This makes calculation of the derivatives simpler 
than for the case of only 2 final state particles. The above equation can be written as 



The Jacobian is now easily evaluated as 


to give 


d(E*,Q*) _ |p/| _ sin#* 
d(Ei,Qi) |p* I sindi 


d 2 a sin#.,* d?a 

dEjidFlji sin#^ dEj^dFlj* ’ 


in agreement with Joachain [19] (p. 30), Dedrick [23], and Byckling [14] (p. 60). This 
expressed purely in terms of 9ji by using 433 to give #j* in terms of Oji. 


5.7 Derivation of d cos 6 JC J d cos 9ji 
Equation (436) gives 


COS 0j c 


—otjc'icil — cos 2 Oji ) ± cos 2 9jiy/D 
7c(! — Pc cos 2 9j l ) 


where 


D = 1 + 7 2 tan 2 9ji( 1 — a 2 c ) . 

We will show that 

d cos 9 j c cos 9ji (ocj C ± \/~D) 2 

dcos9ji 7 2 (1 — (3 2 cos 9 2 l ) 2 (±\/D) 

To simplify the algebra, make the following definitions. 

y = cos 9 jc , 
x = cos 9ji , 

7c = 7 , 
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(462) 

(463) 

(464) 

(465) 

(466) 
may be 

(467) 

(468) 

(469) 

(470) 

(471) 

(472) 





*37. 

Ill 

■37 

(473) 



S 

O 

III 

S 

(474) 



III 

&J 
H |<te 

(475) 

We have 

y 

/ — «7 2 (1 — x 2 ) ± x 2 \J~D 

g 7 2 (1 — /3 2 x 2 ) 

(476) 

where 

/ = 

— cry 2 (l — x 2 ) ± x 2 VD , 

(477) 


3 = 

7 2 (l-/3V) = 7 2 -(7 2 -V, 

(478) 

where we 

have used the result 7 2 /3 

2 = 7 2 — 1. Note that we can also write 



D = 

= 1 + 7 2 (1 — a 2 )(— g — 1) • 

(479) 

Taking derivatives with respect to 

x , we obtain 



D' 

- 

(480) 


r 

= 2c 7 2 x±L^D T 2 * 1 -? 21 , 
L xV(d 

(481) 


±xVDf 

= ± 2 « 7 2 x 2 a/D + 2x 2 D — 7 2 (1 — a 2 ) , 

(482) 


=> ±xVDgf' 

= ± 2 « 7 4 x 2 a/d + 2 x 2 7 2 D — 7 4 (1 — a 2 ) 




=F2a7 4 x 4 >/D — 2x 4 7 2 D + 7 4 (1 — a 2 )x 2 
±2«7 2 x 4 a/D + 2x A D — 7 2 (1 — a 2 )x 2 , 

(483) 


~(±xy/Dg'f ) 

= ±2«7 2 x 2 \fD =F 2a7 2 x 4 \fl) — 2 x a D 




=F 2 a 7 4 x 2 v A D ± 2 a 7 4 x 4 V A D + 2x A 'y 2 D . 

(484) 


Putting this together gives 

±xy/D(gf - g'f) = 2x 2 ^ 2 D + 7 4 (1 - a 2 )(x 2 - 1) 

— 7 2 (1 — a 2 )x 2 ± 2aj 2 x 2 n/D 
= 7 2 x 2 (D ± 2a\fl) + a 2 ) 

= 'y 2 x 2 (a ± VD) 2 , (485) 

finally giving 

gf-g'f _ i 7 2 .T 2 (a ± \/I)) 2 
V g 2 7 4 (1 — /3 2 x 2 ) 2 (±xVD) 
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(486) 


x(a ± \/D ) 2 
7 2 (1 — (3 2 x 2 ) 2 (±\/D) 

Translating back into the original variables gives 

d cos 6j c cos Oji {otj C ± \[D ) 2 

dcosdji 7 2 (1 - /3 2 cos0| z ) 2 (±VD) 

6 Two body final state cross sections 

Three dimensional radiation transport codes require differential cross sections in the lab (space- 
craft) frame. However, calculations are most easily done in the center of momentum or projectile 
frame. The previous section presented general results for Lorentz transforming cross sections, 
energies, angle and momenta from one frame to another. Two body final state cross section 
transformations are much more complicated than three or more bodies, because the energy and 
angle are dependent variables for two body final states. The present section is a culmination 
of this paper, in that explicit expressions are written down for Lorentz transforming 2 - body 
final state cross sections. The key equations (558) - (561) tell one exactly how to transform a 
2 - body final state angular distribution from the center of momentum to the lab (spacecraft) 
frame, with all the relevant kinematic factors written in terms of lab variables. 


6.1 Differential cross sections 

6.1.1 General form of differential cross sections 

The differential cross section is 

S 

da = — |Af| 2 (27r) 4 d<f> 2 (pi + P2]P3,P4) ■ 
Using e = (2tt) 3 2E, the phase space factor is given by 


d®2{pi +P2]P3,P4) 


where it is now understoood that 


X 4 f , \ d 3 P3 d 3 p 4 

<(P.+K-W-M) (2ir)32£ , (2ir) ^ 

s(e 1 + e 2 -e'~ Eijyn 

e 3 e 4 

5(E 1 + E 2 -E , 3 - E 4 )^d\p 4 \dn 4 , 


(488) 


(489) 


P3 = Pi + P2 - P4 , 

so that in the above equation, we have (with | p 4 1 2 = p 2 ) 


E 4 



(490) 


(491) 
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E'z = 


(Pi + P2 - P4) 2 + mj 


(492) 


To eliminate the 5 function, we could now use the Ho-Kim phase space method discussed pre- 
viously [10]. Simplify the above expression by writing 


d|P4 


^P 4 | 

d(m , + E 4 ) 


d(E ' + E 4 ) , 


(493) 


so that d(E 3 + E 4 ) will kill the 5(E\ + E 2 — E 3 — E 4 ) term. It is actually ^ that we will 

evaluate. We use d(E 3 + E 4 ) rather than just (IE 3 because both E ' 3 and E 4 have the term | p 4 1 
present in them. The phase space factor then becomes 


d$2{Pl +P2]P3,P4) 


| P4 1 2 d\p 4 \ 

e 3 e 4 d(E ' 3 + E 4 ) 


d,VL 4 , 


(in general) , 


(494) 


which is identical in form to equation (73), because both phase space factors have a 2 - body 
final state. (The latter equation is obtained from the former via the substitutions 1 — ► 4 and 
2 — ► 3.) The Griffiths phase space method, follows exactly as discussed previously and will give 
the same result. The differential cross section now becomes 


da 

dQ 4 


cS 

64v t 2 F 


\M \ 2 


I P4 1 2 d|p 4 | 

E' 3 E 4 d{E ' 3 + E 4 ) ’ 


(in general) , 


(495) 


which agrees with Ho-Kim [10] (equation 4.57). 


6.1.2 Angular distribution in cm frame 

We now evaluate this expression by evaluating 
differential cross section in the cm frame, where pi + P 2 = 


in a particular frame. Evaluate the 
0 and E 3 = a/ | p 4 1 2 + rri 3 . We obtain 


Also, 


- d(E' 3 + E 4 y 
. d\p 4 \ . 


P 4 


e 3 + e 4 
e’ 3 e 4 


d®2(Pl + P2]P3,P4:)c = 


P4 


d|p 4 | 


e' 3 e4 d(E 3 + E 4 ) 
|p 4 | E' 3 E 4 


d^l 4c 


d^~ti 4 q 


e 3 e 4 E 3 T E 4 
I P4 1 

4(2v r) 6 (E 3 + E 4 ) 


dVt 4c , 


(496) 


(497) 
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which compares to equation (82). The differential cross section is therefore 


da 


S 
4 F 


|Ad| 2 (27r) 4 d<f> 2 (pi +P2]P3,P4) 


S | P4 1 

64v t 2 F E' 3 + E 4 


\M\ 2 dn 4c , 


(498) 


giving 


da 

d£l 4c 


S | P4 1 
64t r 2 F E' 3 + E 4 


\M\ 2 . 


(499) 


Now, use (21) for F with y/s = E 3 + E 4 . The angular distribution in the cm frame is, therefore 
(with |p 4c | = |p 3c | = \pf c \ and |p 2c | = |Pic| = |Pic| in the cm frame), 


da 
d£i 4 C 


S 

647 r 2 s 

cS 

647 r 2 s 


Pfc 

Pic 


\M \ 2 


M2 


where we have used the following results, 


Ip fc 



and 



(500) 


(501) 


The latter expression for is far superior to the former, because the latter is written in terms 
of the single variable s, whereas the former appears to contain two independent variables p ; c 
and p f c . This agrees with references [7] (p. 100), [14] (p. 80), [10] (p. 100), [2] (p. 200). A 
similar formula for is obtained by interchanging 3 4 to give ([14] p. 80) 


da da 
(HI 4 C 


(502) 


The angle is understood to be fli 3 , the angle of particle 3 with respect to particle 1. Written 
explicitly, the above formula becomes 


da da 
dfli 3c dO,24c 


(503) 


When scalar products are encountered such as p 2 • P 4 = |p 4 | |p 2 | cos 0 24 , it becomes necessary to 
specify relative angles such as 0 24 . 


6.1.3 t distribution 

The above results are often written in terms of the Mandelstam invariant variable t defined as 

t = (pi -p 3 ) 2 = (p 2 -P 4) 2 

= - 2 E 2 E 4 + 2p 2 • p 4 + m 2 
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(504) 


giving 


dt = 2|p 2 ||p4|dcos0 2 4 , 


where #24 is the angle between p 2 and P 4 . Because of a minus sign ambiguity to come later, we 
need now to thoroughly establish our convention for a solid angle. Thus, we now present some 
elementary considerations. The solid angle is 

= sin ddOdcp 
= — d cos 9 d(j) 

= — 2nd cos 8 . (506) 

Note that, when integrated this gives the total solid angle. This works as follows. 

[ sin Odd = [— cos0]q = — [cos0]^q = —(—1 — 1) = 2 , (507) 

Jo 

or equivalently, 


sin 8d6 = — d cos 8 = — [cos ( 


IzT 1 = -(- 1 - 1 ) = 2 . 


This gives the total solid angle as 


dQ = 


dip J si nOdQ = 2n J sinOdO 


= —2n 


d cos 8 = An 


Using equations (505) and (506) gives dt = — ^|p 2 ||p 4 |(iU. Using 


da da dQ. c 

dt dQ r dt. 


|P2c||P4c| dQ c ’ 


gives [10] (p. 101), with |p2 C | = |pi c | = |p»c|) and using equation (501) gives 

^ 1 _i A< |2 

dt 647ts |pic| 2 

- 

5 , 


64vrF 2 


|M| 2 . 


Note that all these equations have a minus sign in front, which does not mean the cross section 
is negative, because t takes on negative values. However, most authors, such as de Wit [7] (p. 
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100), Ho-Kim [10] (p. 99), and Byckling [14] (p. 81) do not have a minus sign in such formulas. 
The reason for this is that most authors write 


dU = d cos 9d(j) 

= 27rdcos$, (512) 

instead of equation (506). This is explained as follows. Note that the same answer for equation 
(508) is obtained with 

[ sin 6d6 = + f dcos9 = [cos = (1 1) = 2 . (513) 

Jo J-i 

In other words, we can change the definition of solid angle to 


dU = +27rdcosd, 


(514) 


provided that it is understood that the limits of integration are reversed, as in equation (513). 
This gives the same solid angle as obtained in equation (509) as 


/ 


dil = 


+2-7T / d( cos 6) = 47t 


mi 


(515) 


With this convention, the t integration also changes. Previously, with dVL = — 27rdcos$, the t 
integration would have been fjf dt where to = t(6 = 0) and t n = t{9 = 7 r). However, with our 
new convention, the t integration must be 




+27rdcos# => 



This will be important when we calculate the total cross section, 


a = 



dt . 


(516) 


(517) 


However, because s + t + u = m'f , then u has opposite behaviour to t. Thus, if u integrations 
are being performed, then it must be 


a = 



da 

du 


du . 


With the above considerations, equations (511) then become 


(518) 


da 

dt 


S 


1 


647TS I Pic 


i\M \ 2 
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where by (21) 


cS 

64t tF 2 


\M \ 2 , 


(519) 


A 12 = (s — m 2 — m|) 2 — Am\m\ = 4F 2 . (520) 

The above expressions for ^ are valid in any frame because a and t are invariant. They agree 
with references [1] (equation 34.30) and [14] (p. 81) and [10] (equation 4.62). The invariant 
amplitude often consists of direct (t-channel) and exchange (u-channel) terms, 


M = M d + M e . (521) 

Cross sections, however, always involve the amplitude squared |A4| 2 , and a new interference 
term arises which is the cross term in the square. Thus, 

\M\ 2 = \M d \ 2 + \M e \ 2 + \Mi\ 2 , (522) 


where 


\Mi\ 2 = 2 M d M e . (523) 

Thus, the cross section can be written as a direct, exchange and interference cross section, 

(524) 

-\M e \ 2 , 


day, 

dt 

da e 

dt 

dai 

dt 

The total cross section will be 


5 

167rAl2 

5 


\M d \ 2 


167rAl2 

cS 


cS 


16ttA 12 IMI “ 16ttA 12 2MdMe 


(525) 

(526) 


<7 = 


— Od + <7e + <Xi — 


f to da 

tn ^ 

f to da d 


dt 


f to da. 


*° da.: 


——dt+ / — — dt + / ——dt 


dt 


dt 


dt 


(527) 


When performing integrations to obtain total cross sections, one finds the analysis is much 
simpler if the direct, exchange and interference terms are treated separately. Note that equation 
(3.18) in reference [15] and equation (6) in reference [13] are incorrect; the factor 4 should 
not appear in the denominator. In the notations of those references, the correct equation is 
W = ek\M\ 2 j,. 
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6.1.4 Cross section units 

We consider equation (519) to specifically show how cross section units can be checked. The 
units of the formula (519) for da/dt are verified as follows. On the left hand side, a has units 
of mb or GeV -2 . The Mandelstam variable, t, has units of GeV 2 . Thus, da/dt has units of 
rnb/GeV 2 or GeV -4 . The invariant amplitude will be of the form At = ,_f 2 (for the direct 
term), where t = (p\ — p ^) 2 and g is the coupling constant, and m is the mass of the exchange 
particle. On the right hand side of (519), the coupling g 2 has units of GeV 2 (in scalar theory), 
thus At is dimensionless. A has units of GeV 4 . Thus, the right hand side has units of GeV -4 , 
matching the left hand side. 


6.1.5 Relation between Jy- and ^ 

The differential cross section ^ can be thought of as an angular distribution because t is directly 
related to the scattering angle cos#. Re-writing equation (504) gives 


t. — m 2 — m 2 + 2 E 2 E 4 
2 |P 2 ||P 4 | 


COS #24 • 


(528) 


In the cm frame, Pi + P2 = 0 = P3 + P4 and |pi| = |p2 1 = |pi c | and |p3 1 = | P4 1 = |p/ c |, which 
gives equation (504) as 


t 


(P4 -P 2) 2 

m\ + m| ~ 2^/p 2 C + m 2 \J p 2 c + m| + 2|p/ c ||p ic | cos # 2 4c , 


(529) 


or 


t = ml + ml + — 

A o 


- \J (A12 + 4sm 2 ) (A34 + 4.sm 2 ) + cos #24c 


(530) 


using equation (501). Notice that when mi = m2 = m3 = 7714, this reduces to |p/ c | = |pj c | 
which is used to obtain equation (6.51) in Griffiths [2]. This can also be written [16] (p. 8), [19] 
(p. 31), [14] (pp. 79, 80), 


COS #i 3c = 


s(t — u) + (m 2 — m|)(m. 2 — m 2 ) 
N/A12A34 


(531) 


Byckling [14] (p. 80) and Leon [12] (p. 16) also write equivalent expressions for the lab frame 
angle. Byckling explains how to convert these expressions for the s, t, u channels. Using (514) 
and (530), we obtain 


dt 


y / A12 A34 
2s 


d cos # 


VA12A34 

dil 

47 rs 


(532) 
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4y - is written directly in terms of the angular distribution Using equations (530) and (532), 
we obtain [19] (p. 40) 


da _ \f A12A34 da 
d&24c 47T-S dt ’ 


(533) 


with t given by equation (530), where 624 = O 4 . Since 0 13 = 024 in the cm frame, d 924 /d 9\ 3 = 1, 
therefore 


da da d ^24 da 

dll 24 dfl 13 dll 24 


and 


da _ \/Ai2 A34 da 
dUi 3c 47T-S dt 


(534) 


(535) 


Now, t can also be written as 


t 


(P3 ~Pl) 2 

ml + ml - 2 ^ p 2 c + rafy; pf c + mf + 2|p /c ||p ic | cos0i 3c 


"*1 + m 3 + ^ 


- \/ (A12 + 4smf )(A 34 + 4sru|) + VA12A34 cos 0i 3c 


(536) 


Just as 0 has a minimum and maximum value (i.e. 0 and 7 r), so too does t. The minimum and 
maximum values are denoted 


t 0 = t( 0 c = 0) , 


(537) 


and 


t-K = t( 0 c = 7 r) . 

These are given by [1] (p. 188), [10] (p. 101) 

t 0 {t n ) = (E lc - E 3c ) 2 - (|pic| |p 3c 


(538) 


m{ — ml — m 3 + 777.4 


\/4 s 


- ( Pic T |P3c| 


1 r 
4s - 


{mj -ml~ml + m \) 2 - (\/Ai2 T \^ 34 ^ 


(539) 


where the notation means that to has the — sign and t^ has the + sign, i.e. to = (E\ c — 
E 3 c ) 2 ~ (|pic| - |P3c|) 2 and t n = {E lc - E 3c ) 2 - (|pi c | + |p3c|) 2 - Equation (539) is obtained by 
substituting 0 and ir into equation (536). It is also derived in section 4.4. In plotting an angular 
distribution, cos 9 must not exceed the maximum and minimum values. Similarly, when plotting 
the cross section with t on the horizontal axis, then t must not exceed to or t 
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6.1.6 Lab frame spectral distribution 

In the following derivation, T 4 i is the kinetic energy of particle number 4. The Mandelstam 
invariant t is 


t = (Pi -P 3) 2 = {P2 -PA) 2 , 

= m\ - 2E 4 E 3 + 2pi • p 3 + m\ , 

= m\ - 2E 2 E 4 + 2p 2 • P 4 + ml . (540) 


Previously, we used the 1,3 variables, but now it proves convenient to use the 2,4 variables 
because the lab (target) frame is defined by P 2 = 0 giving E 2 i = m 2 and we obtain 


t = ml - 2m 2 E 4i + , 


or, since E 4 = T 4 + m 4 , 


t = ( m 2 - m 4 ) 2 - 2 m 2 T 4 i , 


which gives 


dt = —2m, 2 dT 4 i . 


(541) 


(542) 


(543) 


de Wit [7] (p. 101) chooses particle number 4 as the target recoil particle. Notice that there 
is no 024 1 appearing in the above two expressions. This makes sense because in the lab frame 
particle 2 is at rest and so it is impossible to have an angle defined relative to particle 2. Using 
equations (542) and (543), we obtain the lab frame spectral distribution of particle 4 ([7] (p. 
101 ) 


do do do 

d'I'.\i = EEki = ~ "' 2 dl ’ 

with t given by equation (542). We also have by (533) 


do 

dT 4 i 


-2 m 2 


do 

dt 


8irm 2 s do 

\/^ 12^34 d$d 4c 
2Trm 2 do 

I Pic I |p/c| d£l 4c 


(544) 


(545) 


The spectral distribution of particle 3 is now easily obtained. Conservation of energy gives 


Eu + E 2 i — E 3 i + E 4l 
= constant. 


(546) 
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Thus. 


dE 3 i = -dE±i , (547) 

dT 3 i = —dT 4 i . (548) 

Alternatively, we could have substituted energy conservation in the lab frame, E 4 = E\ +m 2 — E3 
into equation (542), to obtain 

t. = m\ - ml + 2 m 2 (T 3 ; + m 3 - E lt ) , (549) 


giving 

(550) 

(551) 

with t given by equation (549). 


dt = 2m 2 dT3i 


Thus, 


da da da 

^ = d&u = + "' 2 dl 


6.1.7 Lab frame angular distribution 

To obtain the angular distribution in terms of d 7 , solve (256) for cos$i 3 ; and differentiate, as in 


d(cos 0\3i) 
dE 3 i 


COS 0131 

d{ cos 013 1) 
dT 3l 


E 3 i{Eu + m 2 ) - \{rn\ + ml + 2m 2 Eu + m| - m 2 ) 


lPn||P3/| 


(552) 


E3i{m\ + + ml - m\ + 2 m 2 Eu) - 2m 2 (Eu + m 2 ) 


Eh - m\ {El - m 2 ) 3 / 2 


(lab frame) , 


(553) 


which, upon substitution into equation (551), and with the implied integral over dq i, dTl = 
27t d { cos 0), gives 


da 

m2 

|Pli||P3il 

da 


dQ 3 i 

7 r 

Eu + rn 2 - E31 cos 0 4 31 

dt 



ra 2 

2\J E\i - m\ {El - mi) 3 / 2 

da 


7 r 

E 3 i{ml + ml + ml - m 2 + 2Eum 2 ) - 2ml{Eu + m 2 ) 

dt 


with t given by equation (255). In equation (554), the tall absolute value brackets are necessary. 
E31 is a function of 0 , as shown in equation (265), where in general E31 is a double - valued 
function of 0, and dE3i/d{cos 631) can have both positive and negative values. However, da/did.31 
is positive definite, therefore, must depend only on the absolute value of d(cos 0i 3 i)/ciE 3 /. By 
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substituting for E 3 i from (265), (554) is expressed entirely in terms of #3;. Equation (554) agrees 
with Ho-Kirn [10] (p. 101) , Byckling [14] (p. 80), and deWit [7] (p. 101). Note that if equation 
(553) is substituted into the left hand side of equation (554), a large cancellation of terms on 
the left and right hand sides takes place and one is left back with equation (551), again showing 
the equivalence of (fe and ^ in the lab frame. 

We now present an alternative derivation of the above results. We evaluate equation (495) 
in the lab frame, defined as P2 1 = 0. This means evaluating in the lab frame. E' 3l 

becomes 


E 3 i = y(Pn ~ P<ii ) 2 + ml = sj |pu| 2 + | p 4 z | 2 - 2|pn||p 4 ;| cos 6 + m 3 , 
with 9 = 04 i. This gives 

d(E' 3l + E 41 ) _ . . E' 3 i + E±i{\ — a cos 6 ) 

d|p4zl “ lP4/l i ’ 


a 


P 

|P4i| 


(555) 


(556) 


(557) 


(558) 


(559) 


which agrees with Byckling [14] (p. 80) and Ho-Kirn [10] (equation 4.64). A similar formula for 
is obtained by interchanging 3 4 to give Byckling [14] (p. 80) 


where 


Thus 


da 


S 


P4J 


\Mf 


dQu 647 x 2 F E' u + E 4 ((l — a cos 9) 

Using equation (25) for T) = m2|pu| and E\ + E 2 = Eu + m 2 = E 3 i + E 41 gives 

|M| 2 , 


da 

S 

P4i 

1 

dQ.41 

647r 2 m2 

Pu 

tq 

+ 

to 

1 

pu 

P4i 

E 4 1 COS 641 


da 

S 

P3l 

1 

dQ 3 i 

64:TT 2 m 2 

Pu 

E ll + ™, 2 - 

pu 

P3! 

E 3 1 cos 9 3 i 


I Mf 


(560) 


Upon combining (535) with (554), one obtains 

da 4m2 s 

d^3 1 \J A12 A34 


iPidlPgd 

Eu + m 2 - ^ E 3 i cos 9i 3 i 


4m 2 s 

\A12A34 


E 3 l{m\ + m 2 + m 2 - m\ + 2 Eum 2 ) - 2 m\{Eu + m 2 ) 


da 

dn 13c 

da 

dn 13c 


(561) 
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which provides a way of transforming a cm frame angular distribution to a lab frame distribution. 
This provides an alternative to equation (450). Performing a numerical evaluation of equations 
(450) and (561), which must give the same answer, enables one to perform an excellent test of 
these results. 

7 Conclusions 

This paper has provided a thorough discussion of the transformation of 3 - dimensional differ- 
ential cross sections from the projectile or center of momentum frames to the lab frame. This is 
important because transport codes require all cross sections in the lab frame. Several important 
issues, such as the use of double valued and infinite functions have been discussed extensively. 
Such a detailed and comprehensive treatment does not appear in previous literature. Special 
attention has been paid to writing differential cross sections entirely in terms of lab variables, 
because these are used in 3 - dimensional radiation transport codes. 2 - body final state cross 
sections are much more complicated than 3 - body final states, because the energy and angle are 
dependent variables. Equations (558) - (561) tell one exactly how to transform a 2 - body final 
state angular distribution from the center of momentum frame to the lab (spacecraft) frame, 
with all the relevant kinematic factors written in terms of lab variables. These results will be 
very useful in enhancing space radiation transport codes, such as HZETRN, to be capable of 
fully 3 - dimensional transport. 
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